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Abstract 

Geomagnetic  disturbances  in  the  near  earth  space  environment  can  adversely 
affect  numerous  military  and  Department  of  Defense  (DoD)  systems  and  operations.  To 
improve  the  prediction  accuracy  of  such  disturbances,  the  Global  Assimilation  of 
Ionospheric  Measurements  (GAIM)  working  group  is  spearheading  an  effort  to 
incorporate  near  real-time  ionospheric  measurements  into  the  next  generation  of  space 
environment  forecast  models.  Since  the  model  software  is  designed  to  automate  the  data 
ingest  process,  the  need  arises  to  examine  and  validate  the  quality  of  such  measurements 
before  being  assimilated  into  such  a  model.  One  such  measurement  to  explore,  which  is 
the  focus  of  this  research,  is  the  Defense  Military  Satellite  Program  (DMSP)  measured 
electron  temperature  (Te). 

DMSP  Te  data  were  validated  against  near  simultaneous  incoherent  scatter  radar 
(ISR)  Te  measurements  from  Millstone  Hill,  MA  and  Sondrestrom,  Greenland  for  a  select 
43  conjunctions  between  Winter  1996  and  Summer  2000.  DMSP  Te  measurements  for  a 
given  overpass  were  averaged,  while  ISR  Te  values  were  either  averaged  or  extrapolated, 
depending  on  the  ISR  mode,  to  determine  Te  comparison  values.  In  some  cases, 
instrument  related  anomalies  produced  unreliable  measurements. 

Of  the  37  Millstone  and  six  Sondrestrom  conjunctions  compared,  DMSP  Te 
values  exceeded  ISR  Te  values  by  an  average  of  about  25  percent,  which  is  nearly  three 
times  the  mean  ISR  uncertainty.  Photoelectrons  collected  by  the  DMSP  Electron  Probe 
contaminated  Te  values  particularly  during  solar  minimum.  A  more  comprehensive 
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comparison  extending  to  other  sectors  of  the  DMSP  orbit  is  required  to  determine  if  a 
DMSP  Te  or  ISR  Te  bias  truly  exists. 

Based  on  an  assumed  linear  Te  behavior  at  mid  latitudes,  the  average  DMSP  Te 
random  noise  level  above  Millstone  Hill  was  estimated  at  about  four  percent,  falling  well 
within  the  published  Te  measurement  accuracy  of  10  percent.  This  approach  was 
inappropriate  to  use  at  high  latitudes  due  to  the  high  variability  of  ionospheric 
phenomena  over  short  distances. 
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VALIDATION  AND  ANALYSIS  OF 


DMSP  ELECTRON  TEMPERATURES 
IN  THE  TOPSIDE  IONOSPHERE 


I.  Overview/  Problem  Statement 


Chapter  Overview 

Further  examination  of  Defense  Military  Satellite  Program  (DMSP)  electron 
temperatures  will  immensely  benefit  next  generation  space  weather  forecast  models. 

This  chapter  begins  by  citing  the  motivation  behind  validating  and  assessing  DMSP 
electron  temperature  (Te)  behavior  and  its  impact  on  the  Air  Force  Mission.  Next,  the 
need  to  compare  Incoherent  Scatter  Radar  (ISR)  electron  temperature  measurements  to 
DMSP  electron  temperatures  is  discussed.  The  scope  and  general  approach  of  the 
research  is  then  outlined. 

Introduction 

The  birth  of  the  space  age  in  the  late  1950s  extended  man’s  operational  fi'ontier  to 
include  the  near-earth  space  environment.  With  numerous  assets  routinely  operating  in 
this  environment,  the  Department  of  Defense  (DoD)  relies  upon  accurate  and  timely 
warnings  of  the  implications  of  approaching  geomagnetic  storms  on  its  space-based 
systems  and  on  numerous  current  operations.  After  over  40  years  of  research,  scientists 
have  learned  that  the  near-earth  space  environment’s  composition,  density,  and  dynamics 
can  vary  significantly  with  respect  to  altitude,  latitude,  time  of  day,  season,  solar  cycle. 
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and  geomagnetic  activity.  Although  some  progress  has  been  made,  the  development  of 
space  environment  forecast  models  is  still  in  its  infancy  when  compared  to  its 
tropospheric  coimterpart.  Today’s  space  models,  limited  by  a  delayed  input  of  sparse, 
sporadic  data,  often  act  as  more  of  a  “now”cast,  or  a  given  day’s  ionosphere  “climo”  cast, 
than  as  a  bona-fide  forecast.  Just  as  with  tropospheric  weather,  the  most  reliable  forecast 
models  are  those  created  from  ingesting  a  high  geographical  coverage  of  near  real-time 
measurements  or  observations.  With  a  myriad  of  both  ground  and  space-based  sensors 
coming  on  line,  there  will  soon  be  sufficient  data  coverage  for  assimilation  into  a 
physics-based,  global  ionospheric  model.  Since  the  deployment  of  the  Special  Sensor  for 
Thermal  Ions,  Electrons  and  Scintillations  (SSIES)  instrument  (also  referred  to  as  the 
Topside  Ionospheric  Plasma  Monitor)  on  board  DMSP  spacecraft  in  1987,  the  capability 
has  existed  to  obtain  near  real-time  space  environment  data  in  the  topside  ionosphere  (to 
be  discussed  in  Chapter  II)  at  DMSP  altitude  (approximately  840  km).  The  Global 
Assimilation  of  Ionospheric  Measurements  (GAIM)  working  group  hopes  to  assimilate 
such  real-time  data  along  with  physics-based  algorithms  to  develop  a  reliable  physics- 
based  global  ionosphere  model  within  the  next  five  years  [Schunk  and  Sojka,  19997.  One 
of  the  key  steps  in  developing  such  a  model  is  to  characterize  and  monitor  the  expected 
data  quality,  since  if  bad  data  are  ingested,  the  resulting  ionospheric  specifications  and 
forecasts  will  be  poor.  DMSP  SSIES  measurements  are  one  of  the  several  data  types 
GAIM  plans  to  assimilate;  with  electron  temperatures  being  one  of  the  SSIES  measured 
parameters.  Since  the  model  software  program  is  designed  to  automate  the  data  ingest 
process,  the  need  arises  to  validate  and  investigate  the  quality  of  DMSP  electron 
temperature  data  for  incorporation  into  the  GAIM  algorithms. 
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Impact  on  Air  Force  Mission 

Variations  in  the  near-earth  space  environment  can  adversely  affect  numerous 
military  systems  and  operations.  Increased  geomagnetic  activity  (i.e.  from  solar  flares) 
can  disrupt  high  frequency  radio  communications  and  over-the-horizon  (OTH)  radars, 
alter  Global  Positioning  System  (GPS)  navigation,  and  charge  spacecraft  systems  to  the 
point  of  irreversible  damage.  Due  to  the  ionosphere’s  impact  on  these  systems,  it  is 
critical  that  a  reliable  ionospheric  forecast  model  (i.e.  GAIM)  be  developed.  The  results 
of  this  DMSP  electron  temperature  study  will  directly  benefit  the  GAIM  effort  as  well  as 
Air  Force  Research  Lab’s  (AFRL)  Space  Models  Branch.  A  reliable  GAIM  forecast 
model  could  reestablish  the  standard  for  Air  Force  space  system  protection  and 
optimization. 

Problem  and  Assumptions 

The  focus  of  this  research  involves  investigating  and  statistically  quantifying 
DMSP  electron  temperature  data  for  biases  and  quality  so  that  GAIM  can  adjust  their 
data  ingest  algorithms.  Both  ground-based  and  space-based  sensors  can  measure  such 
ionospheric  parameters  as  electron  and  ion  densities  and  temperatures.  Incoherent  scatter 
radar  (ISR)  is  the  primary  reliable  ground-based  measurement  tool  for  topside 
ionospheric  parameters.  There  are  nine  ISR  sites  worldwide  (Map  of  locations  in  Chapter 
n)  [Millstone  Hill  Website,  2000].  Despite  ISR  data’s  sparse,  intermittent  coverage,  we 
assume  it  to  be  ground  truth  against  which  the  DMSP  data  are  compared.  Incoherent 
scatter  radar  Te  measurements  are  taken  from  the  return  signal  off  of  ionospheric 
electrons.  Since  electron  density  decreases  with  altitude  above  300  km,  ISR  returns  at 
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DMSP  altitude  (840  -  865  km)  are  quite  weak  and  thus  contain  a  low  signal  to  noise  (SN) 
ratio.  Consequently,  a  large  hurdle  to  our  comparison  is  either  missing  ISR  data  or  data 
with  large  error  at  higher  altitudes. 

Meanwhile,  the  SSIES  electron  probe  mounted  on  DMSP  spacecraft  measures  the 
state  of  ionospheric  plasma  at  the  DMSP’s  orbiting  location.  The  electron  probe  records 
continuous  measurements  during  its  101 -minute  sim-synchronous  orbit  around  the  earth. 
Measurements  are  taken  at  all  latitudes,  several  local  times,  all  seasons,  and  all  phases  of 
the  solar  cycle,  but  are  limited  to  DMSP  altitude  (840  -  865  km).  With  SSIES  data 
continually  passed  down  once  per  orbit,  it  is  more  extensive  and  encompassing  than  ISR 
data.  Consequently,  the  challenge  is  to  determine  just  how  dependable,  in  terms  of  its 
validity,  dynamic  range,  and  noise  level,  the  DMSP  electron  temperature  data  really  are. 

Research  Scope  and  General  Approach 

A  previous  draft  document  by  Sultan  and  Rich  [2000]  compares  DMSP  ion 
density  (nj)  measurements  with  near  simultaneous  ISR  ni  measurements  from  Millstone 
Hill  Observatory  in  Massachusetts.  A  recommendation  at  the  end  of  their  paper  was  to 
pursue  a  similar  DMSP  -  ISR  electron  temperature  comparison.  Consequently,  this 
research  involves  validating  and  examining  the  behavior  and  quality  of  DMSP  Te  data 
collected  at  nearly  the  same  time  as  ISR  Te  data  observed  during  the  Plasmaspheric 
Observations  of  Light  Ions  in  the  Topside  and  Exosphere  (POLITE)  campaign  periods 
1-10  that  fell  between  Winter  1996  and  Summer  2000.  SSIES  instruments  mounted  on 
DMSP  spacecraft  F12  -  F15  measure  the  data  covering  this  time  frame.  The  POLITE 
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observation  campaign  is  a  coordinated  effort  by  ISR  sites  worldwide  to  maximize 
simultaneous  coverage  topside  ionospheric  measurements. 

Our  study  concentrates  on  examining  DMSP  electron  temperature  data  collected 
while  orbiting  over  Millstone  Hill,  a  representative  magnetic  mid  latitude  ISR  site,  and 
Sondrestrom,  Greenland,  a  representative  auroral  latitude  ISR  station.  Following  criteria 
similar  to  the  Sultan  and  Rich  nj  study,  coinciding  DMSP  -  ISR  electron  temperature  data 
sets  over  Millstone  Hill  and  Sondrestrom  are  identified  for  comparison  and  validation. 
DMSP  data  within  a  five-degree  circle  of  Millstone  Hill  and  one-degree  circle  of 
Sondrestrom  were  extracted  and  averaged,  while  interpolating  and  extrapolating 
techniques  are  employed  to  determine  a  representative  ISR  electron  temperature  value  for 
comparison.  Due  to  the  complexity  of  Te  behavior  over  high  latitudes  compared  to  mid 
latitudes,  our  efforts  focus  more  on  assessing  the  DMSP  Te  data  above  Millstone  Hill 
than  over  Sondrestrom. 

Once  validated  and  assessed  versus  near  concurrent  ISR  data,  DMSP  Te  data  over 
Millstone  Hill  are  examined  for  random  noise  using  a  simple  linear  fit  approach.  This 
approach  is  not  practical  at  auroral  latitudes  due  to  the  short  scale  lengths  of  ionospheric 
phenomena.  Cases  studies  are  then  presented  which  examine  DMSP  Te  behavior  for 
mid-latitude  ionospheric  phenomena,  potential  instrument  malfunction  or  error,  and 
SSIES  Te  measurement  comparisons  between  spacecraft  with  similar  flyover  times. 

Once  unreliable  cases  are  discarded,  DMSP  Te  trends  are  cited  with  respect  to  noise  level, 
ISR  uncertainty,  solar  cycle,  spacecraft,  and  solar  illumination.  After  these  conclusions 
are  drawn  from  the  results,  recommendations  are  made  for  future  work. 
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II.  Background 


Chapter  Overview 

Accurately  forecasting  variations  in  the  near-earth  space  environment  will  help 
warn  of  potential  impacts  to  various  Air  Force  and  DoD  ground  and  space  based  systems. 
This  chapter  begins  with  some  backgroimd  information  on  the  near-earth  space 
environment,  leading  up  to  electron  temperature  behavior  and  the  topside  ionosphere. 
Events  or  conditions  that  can  alter  topside  electron  temperatures  are  then  discussed. 
These  include  electron  temperature  variations  with  respect  to  solar  illumination,  solar 
cycle,  season,  latitude  and  geomagnetic  activity.  Region  specific  ionospheric  variations 
such  as  sub-auroral  ion  drifts  (SAID),  and  mid-latitude  and  light  ion  troughs  are  then 
briefly  described.  Next,  the  mission  and  characteristics  of  DMSP  spacecraft  are 
introduced,  with  emphasis  on  how  the  Topside  Ionospheric  Plasma  Monitor  (SSIES) 
measures  eleetron  temperatures.  Then  background  details  on  incoherent  scatter  radar 
(ISR)  theory  and  its  electron  temperature  measurement  procedure  are  depicted.  Pros, 
cons,  and  possible  sources  of  error  associated  with  both  DMSP  and  ISR  electron 
temperature  measurement  techniques  are  then  cited.  This  is  followed  by  an  overview  of 
the  purpose  of  the  POLITE  campaign.  The  chapter  eoncludes  with  details  of  some 
previous  DMSP  -  ISR  topside  ionospheric  measurement  comparisons  leading  up  to  this 
DMSP  -  ISR  electron  temperature  validation  study. 
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The  Ionosphere 

The  near-earth  space  environment  contains  a  neutral  and  ionized  atmosphere.  The 
neutral  atmosphere  is  the  lowest,  occupying  a  region  from  the  Earth’s  surface  to  roughly 
60  kilometers  (km).  Just  above  the  neutral  atmosphere,  the  ionosphere  makes  up  the 
lowest  region  of  the  ionized  portion  of  the  earth’s  atmosphere.  The  other  two  ionized 
regions  are  the  plasmasphere  and  magnetosphere.  The  ionosphere  arises  from  the 
interaction  between  solar  radiation  and  the  Earth’s  atmosphere  and  magnetic  field 
[Tascione,  1994]  and  extends  from  about  60  km  to  beyond  1000  km  above  the  earth’s 
surface.  It  consists  of  ions  and  free  electrons  created  primarily  by  photoionization  of 
neutral  atmospheric  species,  mainly  molecular  nitrogen  (N2),  oxygen  (O2),  and  atomic 
oxygen  (O).  Even  though  less  than  1  %  of  the  ionosphere  is  actually  ionized,  there  are 
sufficiently  numerous  free  electrons  and  ions  to  influence  the  propagation  of  radio  waves. 
At  a  given  altitude,  the  ionosphere  reflects  radio  waves  below  the  plasma  frequency  while 
allowing  waves  with  higher  frequencies  to  continue  to  propagate  upward. 

Ionospheric  Density  Structure 

Electron  and  ion  density,  dominant  physical  and  chemical  processes,  and  plasma 
composition,  vertically  distinguish  the  Earth’s  ionosphere.  Electron  density  (ne)  variation 
with  respect  to  altitude  creates  a  basic  layered  structure  at  all  latitudes  [Schunk  and  Nagy, 
2000].  Consequently  the  ionosphere  is  divided  into  four  main  layers  -  D,  E,  Fi  and  F2. 
Figure  1  shows  typical  mid-latitude  daytime  and  nighttime  electron  density  profiles. 
Notice  that  electron  density  exhibits  both  diurnal  and  solar  cycle  variations.  The  four 
distinct  electron  density  regions  develop  due  to  the  differential  solar  absorption 
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Figure  1.  Typical  Mid-latitude  Daytime  and  Nighttime  Electron  Density  Profiles. 
Sxmspot  maximum  and  minimum  are  the  solid  and  dashed  lines,  respectively  [Tascione, 
1994]. 


characteristics  of  the  atmosphere  and  the  resulting  different  physical  and  chemical 
processes  unique  to  each  region.  In  the  D  and  E  regions  (60  -  150  km),  electron  density 
increases  with  altitude  while  chemical  processes  are  the  most  important.  Electron 
densities  continue  to  increase  with  height  in  the  Fi  region  (150-250  km)  as  ion-atom 
interchange  and  transport  processes  start  becoming  important.  The  peak  ion  density  of  up 
to  10^  per  cubic  centimeter  (cm'^)  occurs  in  the  F2  region  at  around  300  km.  Here  the 
ionization  maximum  occurs  due  to  a  balance  between  plasma  transport  and  chemical  loss 
processes  [Schunk  and  Nagy,  2000].  Above  this  F2  peak,  plasma  transport  processes  and 
coulomb  collisions  between  plasma  particles  dominate  relative  to  chemical  processes, 
and  electron  densities  decrease  with  altitude,  eventually  becoming  monotonic.  The 
topside  ionosphere  lies  in  this  region  above  the  F2  peak. 
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The  Topside  Ionosphere 

The  topside  ionosphere  is  the  region  above  the  F2  peak  where  electron  and  ion 
density  starts  decreasing  and  where  remains  the  dominant  ion  [Schunk  and  Nagy, 

2000].  The  topside  ionosphere  ends  and  protonosphere  begins  where  the  lighter  atomic 
ions  (H^  and  He^)  begin  to  outnumber  the  heavier  0"^  ions.  The  primary  reversible  ion- 
atom  interchange  reaction  influencing  O’*"  versus  densities  is: 

0^  +  H  H^  +  0  (1) 

The  dominant  to  light  ion  dominant  transition  altitude  typically  ranges  from 
800  to  1500  km.  Consequently,  the  topside  ionosphere  extends  from  around  300  km  to 
anjwhere  from  800  to  1500  km  (See  Figure  2). 


Figure  2.  Relative  Location  of  the  Topside  Ionosphere  [Schunk  and  Nagy,  2000] 


9 


Historical  Context 


The  existence  of  the  ionosphere  was  established  in  1901  when  Guglielmo 
Marconi  successfully  transmitted  radio  signals  across  the  Atlantic  Ocean  fSchunk  and 
Nagy,  2000].  In  1902,  Arthur  Kennelly  and  Oliver  Heaviside  postulated  that  free 
electrical  charges  in  the  ionosphere  could  reflect  radio  waves  [Ratcliffe,  1967].  The 
following  year,  J.  E.  Taylor  suggested  that  solar  ultraviolet  (UV)  radiation  was  the  source 
of  these  electrical  charges,  implying  solar  control  of  radio  propagation  [Taylor,  1903]. 
The  first  generally  accepted  measurement  techniques  supporting  the  existence  of  the 
ionosphere  occurred  in  1924  by  Appleton  and  Barnett,  and  by  Breit  and  Tuve  ]Schunk 
and  Nagy,  2000].  Breit  and  Tuve’s  pulse  sounding  technique  could  determine  the  height 
of  a  radio  pulse’s  reflection  given  a  signal’s  frequency.  This  remote  sensing  technique 
generates  plots  of  reflection  height  versus  sounder  frequency  called  ionosondes,  which 
are  still  widely  used  today  [Evans,  1975].  Unfortunately  these  early  ground  based 
ionospheric  measurement  techniques  could  only  measure  ionospheric  parameters  up  to 
the  F2  region  peak  at  aroimd  300  km.  Rocket  technology  available  after  World  War  H, 
the  advent  of  the  Space  Age  and  satellite  deployments  in  the  late  1950’s,  and  the 
emergence  of  high  power  pulsed  incoherent  scatter  radar  (ISR)  around  the  same  time 
paved  the  way  for  observing  and  measiiring  physical  characteristics  in  the  topside 
ionosphere. 

Electron  Temperature 

Temperatures  of  microscopic  particles  in  the  ionosphere,  namely  neutrals,  ions, 
and  electrons,  are  kinetic  temperatures,  or  measures  of  energy  of  particles’  thermal 
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motions.  Consequently,  the  average  kinetic  energy  of  particles  in  a  gas  is  directly 
proportional  to  its  temperature.  Temperature  is  a  property  of  a  gas  as  a  whole,  not  of  an 
individual  particle,  thus  the  term  “electron  temperature”  represents  a  temperature 
measurement  of  an  electron  gas.  Thermal  energy  can  be  transferred  throughout  a  gas  via 
conduction,  convection,  or  radiation.  Thermal  conduction,  or  the  transfer  of  thermal 
energy  by  contact,  is  the  dominant  means  of  thermal  energy  displacement  in  the  topside 
ionosphere.  Photoelectrons  from  solar  extreme  ultraviolet  (EUV)  radiation  provide  the 
main  source  of  energy  for  thermal  electrons  at  all  latitudes.  These  photoelectrons  directly 
transfer  energy  to  topside  thermal  electrons  via  coulomb  collisions,  thus  heating  the 
electrons.  Higher  electron  temperatures  can  generally  be  expected  at  high  latitudes, 
where  precipitating  auroral  electrons  provide  an  additional  source  of  energy.  Topside 
thermal  electrons  lose  energy  and  thus  cool  via  coulomb  collisions  with  ions  and 
downward  thermal  conduction  [Schunk  and  Nagy,  2000] . 

Electron  Temperature  vs.  Electron  Density  Relationship 

For  a  given  electron  heating  rate,  the  electron  temperature  is  inversely 
proportional  to  the  electron  density.  Lower  electron  densities  will  lead  to  greater  electron 
temperatures  due  to  the  greater  thermal  energy  available  per  particle,  and  less  coulomb 
coupling  with  ions  [Moffett  and  Quegan,  1982].  Conversely,  a  higher  concentration  of 
electrons  equates  to  more  electrons  sharing  a  given  amount  of  thermal  energy,  thus  the 
equivalent  electron  kinetic  temperature  is  lower.  This  inverse  Te  -  iie  relationship  does 
not  always  hold,  as  the  amount  of  thermal  energy  does  vary  with  respect  to  solar  angle 
(Discussed  on  page  14),  solar  cycle,  and  season.  Also  the  plasma  recombination  rate. 
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which  reduces  Ue,  depends  non-linearly  on  Te,  thus  a  slight  increase  in  Te  can  trigger  a 
much  larger  reduction  in  ne  [Sojka,  private  correspondence,  8  Feb  01]. 

Topside  Te  Profile  -  Diurnal  Variation 

The  main  source  of  ionization  and  thermal  energy  for  the  ionosphere  is 
photoionization  [Schunk  and  Nagy,  2000].  Thus  when  the  ionosphere  is  sunlit  and  as  the 
solar  zenith  angle  (see  page  14)  decreases,  ionization  producing  photoelectrons,  and 
photon  thermal  energy  (hv)  increase.  Typical  mid  latitude  Te  values  in  the  topside 
ionosphere  range  from  1200  K  (.10  eV)  at  night  to  3500  K  (.30  eV)  during  the  day 
[Sultan  and  Rich,  Draft  Document,  2000].  This  Te  diurnal  variation  is  chiefly  due  to  the 
photoelectron  heat  source  being  shut  off  at  night.  At  sunrise,  Te  and  ne  increase  rapidly 
as  a  result  of  photoelectron  heating  and  photoionization  [Schunk  and  Nagy,  2000].  Te 
usually  levels  off  towards  midday,  and  then  decreases,  as  does  rie,  towards  sunset  as  the 
photoelectron  heat  source  disappears.  Te  values  at  topside  altitudes  normally  exceed  the 
temperatures  of  neutrals  (Tn)  and  ions  (Ti).  This  is  primarily  based  on  a  conservation  of 
momentum  and  mass  argiunent,  with  electrons  having  the  smallest  mass,  and  thus 
exhibiting  the  fastest  thermal  motion  and  highest  temperature  of  the  three  species.  The 
higher  Te  values  can  also  be  attributed  to  the  fact  that  ion-electron  and  electron-neutral 
coupling  are  low  compared  to  ion-neutral  coupling.  Figure  3  depicts  characteristic  day 
and  nighttime  Tn,  Ti,  and  Te  profiles  up  to  800  km  over  Millstone  Hill  Observatory  while 
Figure  4  represents  an  average  mid  latitude  temperature  profile  up  to  4000  km.  Notice 
how  Ti  comes  closer  to  but  never  exceeds  Te  as  altitude  increases. 


12 


Figure  3.  Typical  Te,  T,,  and  Profiles  for  Daytime  (left)  and  Nighttime  (right) 
Ionospheres  over  Millstone  Hill  on  23-24  March,  1970  [Schunk  and  Nagy,  2000] . 


Figure  4.  Temperature  Profile  of  a  Typical  Daytime  Mid-latitude 
Ionosphere.  Altitude  (km)  represents  the  vertical  axis. 
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Topside  Te  Variations 

The  Te  profile  above  400  km,  where  thermal  conduction  dominates,  usually 
increases  with  altitude  in  response  to  a  downward  heat  flow  from  the  magnetosphere 
[Schunk  and  Nagy,  2000].  Topside  electron  temperatures  typically  undergo  solar  zenith 
angle  (SZA)  and  seasonal  variations.  The  ionosphere’s  seasonal  variation  is  mainly 
related  to  a  solar  zenith  angle  change.  The  SZA  is  the  angle  away  fi'om  which  the  sun’s 
rays  are  directly  overhead.  The  solar  zenith  angle  of  a  satellite  orbiting  the  Earth’s 
ionosphere  can  be  approximated  by  the  angle  formed  between  the  sun  -  center  of  the 
Earth  line  and  the  center  of  the  earth  -  satellite  line  (See  Figure  5). 


Based  on  an  average  DMSP  orbiting  altitude  of  850  km,  trigonometry  was  used  to 
calculate  the  terminator  solar  zenith  angle  for  DMSP  spacecraft  to  be  1 1 8  degrees.  A 
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zero  degree  zenith  angle  implies  a  directly  overhead  sun,  thus  a  smaller  solar  zenith  angle 
means  more  direct  sunlight  and  more  frequent  photoionization  at  a  given  altitude.  Thus 
for  a  fixed  amount  of  energy,  a  smaller  solar  zenith  angle  results  in  increased  iie  and  a 
corresponding  Te  decrease.  Conversely  a  larger  zenith  angle,  as  found  in  the  winter 
hemisphere,  suggests  less  photoionization  at  a  given  altitude,  lower  electron  densities  and 
higher  electron  temperatures  for  a  fixed  amount  of  energy.  The  overall  Te  -  rie 
correlation  is  not  exclusively  tied  to  solar  zenith  angle  variation.  Neutral  circulation 
patterns  and  thermospheric  coupling  with  the  ionosphere  create  a  seasonal  anomaly, 
whose  effects  can  outweigh  SZA  influence  on  Te  and  iie  [Schunk  and  Nagy,  2000]. 

Electron  temperatures  also  vary  with  respect  to  solar  cycle.  The  solar  cycle  is  an 
approximate  1 1-year  variation  in  the  amount  of  sunspot,  solar  flare,  and  geomagnetic 
storm  activity.  Solar  maximum  corresponds  to  maximum  solar  activity,  while  solar 
minimum  signifies  the  low  point  of  solar  activity.  At  solar  maximum,  the  solar  extreme 
ultra-violet  (EUV)  radiation  fluxes  are  largest,  leading  to  increased  photoionization 
resulting  in  higher  electron  densities  and  lower  electron  temperatures.  Conversely 
electron  temperatures  are  higher  at  solar  minimum  due  to  lower  electron  densities.  The 
solar  10.7  cm  radio  flux  acts  as  a  reasonable  estimate  for  the  solar  EUV  flux  [Rich, 
DMSP  Website,  2000],  Thus  tracking  and  plotting  the  10.7  cm  flux  trend  over  several 
years  can  reveal  the  solar  cycle  pattern  (See  Figure  14  in  Chapter  m). 

Topside  electron  temperatures  also  usually  vary  with  geomagnetic  latitude.  As 
mentioned  earlier,  photoelectrons  at  topside  altitudes  transfer  energy  to  thermal  electrons 
by  thermal  conduction  as  they  stream  along  magnetic  field  lines.  Photoelectrons  have 
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more  of  an  opportunity  to  transfer  energy  to  thermal  electrons  along  longer,  higher 
latitude  field  lines,  than  along  shorter,  lower  latitude  ones  [Schunk  and  Nagy,  2000], 

Overall,  electron  temperatures  at  high  latitudes  are  characteristically  subject  to  the 
most  variations.  At  high  latitudes,  strong  electric  fields  arising  from  the  interaction  of  the 
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Earth’s  magnetic  field  with  the  interplanetary  magnetic  field  (IMF)  produce  swift  Ex — j 
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ion  and  electron  drifts  [Rishbeth,  1988],  These  drifts  increase  the  likelihood  of  plasma- 
neutral  collisions  that,  in  turn,  induce  frictional  heating  and  can  elevate  ion  and  electron 
temperatures.  In  addition,  high  latitude  geomagnetic  field  lines  can  connect  to  the 
Earth’s  magnetotail,  allowing  for  particle  precipitation  from  the  tail  to  enter  the  topside 
ionosphere.  This  incoming  auroral  precipitation  can  act  as  a  source  of  bulk  heating  and 
ionization.  Electron  temperatures  increase  as  a  result  of  these  precipitating  electrons 
transferring  energy  to  existing  thermal  electrons  via  coulomb  collisions  that  in  turn  create 
energetic  secondary  electrons  through  ionization  [Schunk  and  Nagy,  2000] .  On  the  flip 
side,  the  high  latitude  “open”  field  lines  also  allow  ions  and  electrons  to  escape  from  the 
topside  ionosphere  as  well.  These  aforementioned  processes  at  high  latitudes  all  produce 
short  scale  length  variations  on  the  topside  ionosphere’s  electron  temperature  and  density 
structure. 

Mid  Latitude  Ionospheric  Phenomena 

At  geomagnetic  mid-latitude  regions,  where  the  DMSP  -  ISR  electron 
temperature  comparison  above  Millstone  Hill  takes  place,  Te  behavior  is  normally  quite 
smooth  compared  with  high  latitude  Tg  variations,  but  still  can  deviate  due  to  distinctive 
ionospheric  phenomena.  Some  imiquely  mid-latitude  ionospheric  phenomena  that  can 


16 


alter  mid  latitude  topside  electron  temperatures  are  sub-auroral  ion  drifts  (SAID),  mid 
latitude  troughs  and  light  ion  troughs. 

A  SAID  event  is  a  latitudinally  narrow  region  (from  0.1-2  degrees)  of  a  rapid 
westward  ion  drift,  located  in  the  evening  sector  on  the  equatorward  edge  of  the  auroral 
oval  [Anderson  et  al,  1993].  These  ion  drifts  can  reach  as  high  as  4  km  s‘\  A 
corresponding  Ti  spike  and  ne  drop  also  occur  during  a  SAID.  Although  SAIDs  are  most 
prominent  at  F2  altitudes  (~  400km),  upward  thermal  conduction  could  result  in  elevated 
Te  values  at  DMSP  altitude.  SADDs  typically  occur  during  the  recovery  phase  of  a 
geomagnetic  substorm  and  last  less  than  three  hours  [Anderson  et  al,  1993].  During  a 
SAID  event,  large  poleward-directed  electric  fields  develop  that  can  produce  westward 
ion  drifts.  These  enhanced  electric  fields  drive  ions  to  collide  with  neutrals  and  then 
recombine  with  electrons.  This  chain  of  reactions  depletes  total  ion  and  electron  densities 
while  at  the  same  time  increasing  ion  and  electron  temperatures.  The  most  notable  duo  of 
reactions  associated  with  this  process  is: 


0^  +  N2^N  +  N0^ 

(2) 

NO^  +  e'  ^  N  +  0 

(3) 

Due  to  the  large  thermal  conductivity  of  the  topside  ionosphere,  fiictional  heating 
associated  with  these  ion  -  neutral  collisions  can  contribute  to  elevated  electron 
temperatures  at  DMSP  altitudes. 

Turning  to  the  next  mid-latitude  feature  that  can  affect  electron  temperature, 
according  to  Schunk  et  al.  [1976],  the  mid  latitude  trough  (MLT)  is  a  region  of  low 
electron  density  typically  found  just  equatorward  of  the  night  side  auroral  oval.  Due  to 
the  inverse  relationship  between  Te  and  ne  for  a  given  amount  of  energy,  a  corresponding 
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region  of  elevated  Te  may  be  expected  [Moffett  and  Quegan,  J982J.  The  MLT  is  most 
pronounced  during  fall  and  winter  in  the  nighttime  ionosphere,  especially  near  magnetic 
midnight.  Looking  at  characteristic  MLT  morphology,  the  base  of  the  trough  is  about 
four  to  five  degrees  of  latitude  wide,  with  the  poleward  trough  wall  being  steeper  than  the 
equatorward  wall  [Schunk  et  al.,  1976].  Electron  densities  in  the  base  on  the  trough  can 
be  on  the  order  of  five  to  ten  times  lower  than  those  on  either  side  of  the  trough. 

Although  the  MLT  is  most  evident  at  lower  F-region  altitudes,  its  characteristics 
are  still  noticeable  at  topside  altitudes.  The  processes  that  can  deplete  electron  densities 
and  thus  contribute  to  the  formation  of  the  MLT  include  the  absence  of  photoionization  at 
night,  plasma  escape  via  the  polar  wind,  plasma  convection  and  ion  chemistry  [Schunk  et 
al,  1976].  Both  electrons  and  are  removed  from  the  topside  ionosphere  by  way  of  the 
ion  chemistry  equations  (2)  and  (3)  on  the  previous  page,  thus  supporting  MLT  ne  trough 
development  [Schunk  et  al,  1976]:  The  impact  that  a  particular  process  has  on  MLT 
trough  formation  depends  upon  a  combination  of  factors  such  as  local  time,  season,  the 
size  of  the  auroral  oval,  and  the  level  of  past  and  present  magnetospheric  disturbances 
[Schunk  et  al,  1976]. 

Another  mid-latitude  ionospheric  featme  that  can  be  observed  at  topside  altitudes 
is  the  light  ion  trough  (LIT).  The  LIT  is  characterized  by  a  pronounced  decrease  of  the 
light  ion  (H^  and  He"^)  densities,  sometimes  with  little  or  no  corresponding  decrease  in  the 
electron  density  (rie).  and  He"^  densities  may  drop  by  as  much  as  two  orders  of 
magnitude  within  5-10  degrees  of  latitude  [Taylor,  1972];  see  Figure  6.  Although  the 
LIT  can  persist  during  day  and  night,  it  is  most  pronounced  during  all  seasons  at  night, 
and  during  da)dime  in  the  winter  hemisphere.  The  steepest  gradient  heading  into  the  base 
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of  the  LIT  is  typically  located  near  60  degrees  geomagnetic  latitude,  with  a  latitudinal 
extent  of  5  -  10  degrees,  but  can  move  equatorward  and  deepen  in  response  to 
geomagnetic  storms  [Taylor,  1972].  Thus  Millstone  Hill  Observatory,  located  at  53.2 
degrees  geomagnetic  north  latitude,  could  occasionally  be  within  the  equatorward  edge  of 
a  LIT.  Taylor  (1972)  adds  that  the  base  of  the  LIT  appears  to  mark  the  average  position 
of  the  plasmapause,  or  the  boundary  separating  the  more  dense  plasma  (H"^)  corotating 
with  the  Earth  from  the  less  dense  plasma  affected  by  magnetospheric  electric  fields.  At 
the  base  the  LIT,  and  He^  densities  can  drop  by  as  much  as  two  orders  of  magnitude 
to  as  low  as  10  ions  /  cm  [Taylor,  1972].  No  firm  conclusion  was  drawn  showing  any 
sort  of  Te  fluctuation  in  response  to  a  light  ion  trough.  If  the  light  ions  (H"^  or  He"^)  are 
the  dominant  ions  at  the  onset  of  a  LIT,  electron  temperature  can  increase  in  response  to 
LIT  development.  On  the  other  hand,  Te  would  show  very  little  variation  in  response  to  a 
LIT  if  the  light  ions  are  the  minor  ions. 


Figure  6.  A  LIT  Measured  by  the  OGO-2  Satellite.  Light  ion  densities  can  drop 
an  order  of  magnitude  with  a  few  degrees  latitude  [Shunk  and  Nagy,  2000], 
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Incoherent  Scatter  Radar 


Incoherent  scatter  radar  (ISR)  is  one  of  the  most  powerful  ground  based  remote 
sensing  techniques  able  to  provide  consistent  measurements  of  the  topside  ionosphere. 

An  ISR  sites  transmits  shorts  bursts,  or  pulses,  of  electromagnetic  energy  to  determine 
the  incoherent  scatter  echo  by  electrons  in  the  ionospheric  plasma.  The  length  of  the 
pulse  determines  the  altitude  resolution  of  the  return.  The  ISR  principle  is  based  upon  J. 
J.  Thomson’s  1906  discovery  showing  that  single  free  electrons  are  capable  of  weakly 
scattering  electromagnetic  waves  [Schunk  and  Nagy,  2000].  Thus  incoherent  scatter  is 
often  called  “Thomson  scatter”  on  his  behalf  The  term  “incoherent  scatter”  refers  to  the 
scattered  radar  return  signal  from  ionospheric  electrons  exhibiting  random  thermal 
motions.  Since  the  electrons  move  at  varying  thermal  velocities,  the  incoherent  scatter 
return  echo  portrays  a  Doppler  broadened  range,  or  spectrum,  of  frequencies  surrounding 
the  transmitter  frequency  [Millstone  Hill  ISR  Tutorial,  2000].  An  extremely  powerful 
radar  system  is  required  to  detect  the  weak  incoherent  scatter  echo.  This  return  echo  can 
be  as  much  as  18  orders  of  magnitude  weaker  than  the  outgoing  signal  (picowatts  vs. 
megawatts). 

Determining  ISR  Electron  Temperature 

Despite  being  much  weaker  than  the  original  signal,  the  ISR  return  signal  contains 
valuable  information  about  ion  and  electron  temperature,  composition,  and  velocity. 
Determining  a  reliable  electron  temperature  from  the  weak  incoherent  scatter  signal 
involves  making  some  key  assumptions.  First,  it  is  assiuned  that  the  ionospheric  plasma 
exhibits  a  Maxwellian  velocity  distribution,  where  particle  collisions  dominate,  and  ions 
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and  electrons  move  at  the  mean  thermal  speed  (( v,. )  -  see  Equation  6)  [Schunk  and  Nagy, 
2000].  Secondly,  although  incoherent  scatter  comes  from  electrons,  ions  influence 
electron  motions  as  long  as  the  radar  wavelength  (k)  is  much  greater  than  the  Debye 
wavelength  [Hargreaves,  1992],  Schunk  and  Nagy  [2000]  define  the  Debye  length  as 
the  minimum  distance  over  which  a  plasma  can  exhibit  collective  behavior.  The  Debye 
length  is  defined  as: 


Xd= 
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where 


So  =  Permittivity  of  free  space 
k  =  Boltzmann  constant 
Te  =  Electron  temperature 
Me  =  Electron  density 
e  -  Elementary  charge 

The  incoherent  scatter  spectrum  includes  two  weak,  narrow  “electron  lines,”  and  a 
broader,  more  vivid  double  peaked  ion  line  power  spectrum  [EISCAT  Website,  2000]. 
Most  ISR  sites  are  designed  primarily  to  observe  and  analyze  the  ion  line  spectrum  to 
determine  electron  temperature.  This  is  based  on  the  assumption  that  k  »  Xd,  ions 
influence  electron  motions,  implying  that  the  major  portion  of  the  ISR  return  is 
concentrated  in  the  ion  line  [Schunk  and  Nagy,  2000].  Given  the  Debye  length  definition 
above,  the  Debye  length  increases  as  electron  densities  decrease.  As  the  Debye  length 
approaches  or  exceeds  the  radar  wavelength,  scattering  comes  from  individual  electrons, 
resulting  in  a  Gaussian  shaped  electron  line  return  instead  of  the  ion  line.  This  could 
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disrupt  the  assumed  technique  of  determining  Te  from  the  dual  humped  ion  line  return. 
The  twin  ion  line  peaks  correspond  to  Doppler  shifts  (Afi)  due  to  ions  approaching  the 

radar  at  the  mean  thermal  speed  (v,. ).  Figure  7  below  shows  such  a  typical  ion  line 
backscatter  spectrum  where  Te  =  Tj. 


T,  .  T,  Afi  »  ^ 


Figure  7.  Typical  Doppler  Broadened  Backscattered  Ion  Line  Power  ISR  Spectrum  in  the 
Case  Where  Te  =  Ti  [Evans,  1975]. 


The  ion  line  Doppler  shift,  Af,  dependent  upon  v,. ,  can  be  expressed  as  (in  Hz): 


where 
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(6) 
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Substituting  for  the  mean  ion  thermal  speed  ( v,. )  )delds: 
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where 


k  =  Boltzmann  constant 
Ti  =  Ion  temperature 
Mi  =  Ion  mass 
X  =  Radar  wavelength 

Since  ion  mass  is  part  of  the  Doppler  shift  expression,  the  ionospheric  ion  composition  is 
instrumental  in  determining  A/i,  Tj,  and  ultimately  Te.  Given  the  above  relationships,  ion 
temperature  is  determined  from  the  width  of  the  twin  peaks  (as  displayed  in  Figure  7), 


while  the 


ratio  is  extracted  from  the  ion  line  spectrum  shape  [Evans,  1975].  A 


wider  ion  line  spectrum  indicates  higher  ion  temperatures,  and  a  larger  Doppler 
frequency  shift  leads  to  a  more  shallow  ion  line  spectrum  in  the  center  with  more 
pronoimced  peaks  toward  the  edges  [Millstone  Hill  ISR  Tutorial,  2000].  Figure  8  shows 


how  the 


ratio  affects  the  positive  Doppler  shift  wing  of  the  incoherent  scatter  echo 


ion  line  power  spectrum.  Once  the  ion  temperature  and 
electron  temperature  can  be  easily  calculated. 


ratio  are  identified,  the 
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Figure  8.  Positive  Doppler  Shifted  ISR  Power  Spectra  for  1!JT\  =  1  to  4 
[Hargreaves,  1992] 


Incoherent  Scatter  Radar  Sites 

Electron  temperature  data  from  two  of  the  nine  ISR  sites  worldwide  were  used  for 
this  research  -  Millstone  and  Sondrestrom.  Millstone  Hill  observatory  in  Massachusetts 
served  as  a  representative  mid-latitude  site,  while  Sondrestrom,  Greenland  observatory 
acted  as  a  symbolic  high-latitude  site.  Figure  9  shows  a  map  of  all  nine  locations,  and 
Table  1  reveals  coordinate  and  transmission  fi-equency  details  for  these  two  ISR  sites. 
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Figure  9.  Worldwide  ISR  Locations  [Millstone  Hill  Website,  2000] 


Table  1.  ISR  Sites  Used  in  this  Study  [Schunk  and  Nagy,  2000] 


ISR  Site 

Geographic 

Latitude 

Geographic 

Longitude 

Geomagnetic 

Latitude 

Elevation 

Angle 

Transmission 

Frequency 

(MHz) 

Millstone 

Hill,  MA 

42.6°  N 

71.5°  W 

53.2°  N 

O 

00 

00 

440 

Sondrestrom, 

Greenland 

67.0°  N 

51.0°  W 

73.1°  N 

o 

o 

00 

1290 

Sondrestrom 
“Up  B” 

(@  850km) 

65.9°  N 

48.8°  W 

N/A 

00 

o 

o 

1290 

The  Atmospheric  Sciences  Group  (ASG)  of  the  Massachusetts  Institute  of 
Technology’s  Haystack  Observatory  operates  the  Millstone  Hill  ISR  observatory 
[Milstone  Hill  Website,  2000].  ISR  measurements  began  at  Millstone  Hill  in  I960.  The 
ISR  instruments  there  consist  of  two  2.5  MW  440  MHz  transmitters.  One  is  a  fully 
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steerable  46-meter  antenna,  while  the  other  is  a  zenith  directed  68  meter  fixed  antenna 


(see  Figure  10). 


Figure  10.  Millstone  Hill  Fixed  Dish  (left)  and  Steerable  Dish  (right)  ISR  Antennas 
[Millstone  Hill  website,  2000] 

Millstone  Hill’s  location  makes  it  ideal  to  study  mid-latitude  (sub-auroral)  ionospheric 
phenomena  and,  on  occasion,  the  aurora.  Typical  Millstone  Hill  ISR  pulse  length 
measurements  are  at  410, 640, 1000,  and  2000  microseconds  (ps). 

SRI  Intemational,  based  in  Menlo  Park,  CA,  operates  the  Sondrestrom  ISR  station 
in  cooperation  with  the  Danish  Meteorological  Institute.  The  site  has  been  in  operation 
since  1983.  ISR  data  are  collected  with  a  32-meter,  fully  steerable  dish  antenna  with  a 
3.5  MW  1290  MHz  transmitter  (see  Figure  1 1).  Sondrestrom’s  location  is  ideal  for 
measuring  the  auroral  and  sub-auroral  ionosphere. 
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Figure  1 1 .  Sondrestrom’s  32  m  Steerable  Dish  [Sondrestrom  website,  2000] 


The  DMSP  Program 

Defense  Meteorological  Satellite  Program  (DMSP)  satellites  are  a  series  of  sun 
synchronous  polar  orbiting  satellites  that  orbit  the  Earth  at  altitudes  ranging  from  840  - 
865  km.  Each  satellite  has  an  orbital  period  of  approximately  101  minutes.  A  minimum 
of  two  spacecraft  is  usually  in  orbit  simultaneously.  The  spacecraft  are  in  orbits  fixed  in 
local  time  with  equatorial  crossing  times  near  0600  (descending)  and  1800  (ascending), 
and  near  0900  (descending)  to  2100  (ascending).  Table  2  shows  equatorial  crossing 
times  for  the  four  spacecraft  used  in  this  research  (F12  -  F15). 

The  primary  mission  of  DMSP  spacecraft  is  to  observe  tropospheric  weather.  The 
secondary  mission  of  DMSP  is  to  monitor  the  near-Earth  space  environment.  Since 
1987,  a  Topside  Ionospheric  Plasma  Monitor  consisting  of  Special  Sensors  for  Thermal 
Ions,  Electrons,  and  Scintillation  (SSIES)  has  been  mounted  on  DMSP  spacecraft  to 
measure  thermal  plasma  parameters  along  the  satellite  flight  path. 
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Table  2.  Recent  DMSP  Equatorial  Crossing  Local  Times  (Ascending  /  Descending) 
[Rich.  DMSP  Website.  2000] 

Satellite  Nominal  Local  Time  of  Node  Actual  Local  Time  of  Node 


In  particular,  the  SSIES  package  measures  ion  and  electron  density,  temperature,  sensor 
potential,  ion  drift  velocity,  and  He^,  and  ion  composition.  Figure  12  below  is  a 
schematic  of  a  DMSP  spacecraft  in  its  deployed  configuration. 


Figure  12.  SSBES  Sensors  Mounted  on  a  DMSP  Spacecraft  [Rich.  1994] 
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The  Langmuir  Electron  Probe 

The  specific  SSIES  instrument  that  measures  electron  temperature  is  the  electron 
probe  (EP).  This  is  typically  a  spherical  or  cylindrical  Langmuir  probe  engineered  to 
collect  charged  particles  in  surrounding  environmental  plasma.  The  EP  mounted  on  the 
DMSP  spacecraft  is  a  spherical  Langmuir  probe  designed  to  measure  anticipated  electron 
temperatures  at  DMSP  altitude  at  all  latitudes,  seasons,  and  phases  of  the  solar  cycle. 
With  this  in  mind,  the  EP  measure  electrons  in  the  temperature  range  of  500  -  9000  K  to 
within  10  percent  accuracy  [Rich,  1994].  The  EP  is  mounted  on  a  2.5-meter  boom 
extending  out  from  the  body  of  the  spacecraft  to  isolate  its  potential  fi'om  that  of  the 
spacecraft  and  to  minimize  collection  of  unwanted  photoelectrons  from  the  body  of  the 
spacecraft  (see  Figure  13).  The  electron  probe  weighs  0.2  lbs  and  consists  of  an  outer 
gridded  sphere  of  diameter  2.25  inches,  and  a  solid,  inner  sphere  of  diameter  1.75  inches. 


Figure  13.  Electron  Sensor  Mounted  on  the  End  of  its  Boom  [Rich,  1994] 
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There  is  a  slight  variation  of  the  SSIES  electron  probe  design  compared  to  a 
standard  Langmuir  probe.  A  standard  Langmuir  probe  consists  of  a  standard  spherical 
conducting  surface  that  is  directly  exposed  to  the  plasma  and  collects  ions  and  electrons 
as  the  probe  potential  is  varied  [Rich,  1994].  Meanwhile,  the  SSIES  electron  probe  has 
an  additional  outer  spherical  grid  that  acts  as  an  electrostatic  screen  filtering  out  ambient 
ions  from  reaching  the  charged  particle  collecting  inner  sphere.  The  inner  sphere’s 
voltage  is  electrostatically  biased  by  +20  V  to  create  a  positive  potential  (Vbias)  with 
respect  to  the  outer  grid  and  spacecraft  body  to  attract  electrons  through  the  grid  and 
repel  thermal  ions.  Rich’s  SSIES  Users  Guide  [1994]  states  that  although  the  outer  grid 
repels  the  thermal  ions,  other  imwanted  charge  carriers  such  as  auroral  electrons  and 
photoelectrons  ejected  from  the  wires  of  the  grid  are  still  collected,  possibly 
contaminating  Tg  measurements  (see  “Sources  of  Error”  section,  pg.  32).  According  to 
Fred  Rich’s  DMSP  web  site  [2000],  the  EP  completes  one  measurement  cycle  every  8 
seconds,  consisting  of  the  potential  on  the  EP  grid  being  swept  from  -4  V  to  +4  V  and 
vice  versa  with  respect  to  Vbias-  Each  4  seconds  of  data  produces  one  measurement  of 
electron  temperahrre  and  density.  The  variation  of  collected  current  versus  applied 
voltage  helps  determine  electron  temperatures  (see  next  section). 

Electron  Temperature  Determination 

Electron  temperature  is  determined  from  evaluating  the  logarithmic  slope  of  a 
best  fit  line  to  the  electron  retarding  region  of  a  curve  plotting  EP  collected  electron 
current  (Ig)  versus  EP  potential  relative  to  the  spacecraft  potential  [Schunk  and  Nagy, 
2000].  Varying  the  EP  voltage  and  measuring  the  resulting  EP  eurrent  generates  this  EP 
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current  versus  EP  potential  curve.  The  electron  retarding  region  is  where  the  EP  potential 
relative  to  the  spacecraft,  or  Vp,  is  less  than  the  plasma  potential  relative  to  the  spacecraft, 
Vo.  Assuming  a  Maxwellian  distribution  of  the  electron  gas,  the  electron  current 
collected  by  the  EP  in  the  electron  retarding  region  is  given  by  [Schunk  and  Nagy,  2000] : 


where 
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le  =  EP  Electron  Cmrent 
e  =  Electron  Electric  Charge 
He  =  Electron  Density 
A  =  Area  of  Gridded  Sphere 
k  =  Bolt2mann  Constant 
Te  =  Electron  Temperature 
/We  =  Electron  Mass 

Vp  =  Probe  Potential  Relative  to  Spacecraft  Potential 
Vo  =  Plasma  Potential  Relative  to  Spacecraft  Potential 


(8) 


Taking  the  logarithm  of  (8)  gives: 
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where 
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Taking  the  derivative  of  (9)  with  respect  to  the  EP  potential  gives: 


d(\ogIe)  e 


dVp 


kTe 
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After  eliminating  data  points  caused  by  photoelectrons  and  hot  electrons,  a  linear  fit  is 
applied  to  the  electron  retarding  portion  of  the  log  le  vs.  Vp  curve.  The  left  hand  side  of 
(10)  represents  the  slope  of  this  line.  Rearranging  terms  and  solving  for  Te  (in  Kelvin) 
yields: 

£ 

k  d(\ogIe) 


Pros,  Cons,  and  Sources  of  Error 

There  are  several  differences,  advantages  and  disadvantages  to  each  method  of 
measuring  parameters  in  the  topside  ionosphere.  A  major  difference  between  the  DMSP 
Te  and  ISR  Te  measurement  techniques  is  that  the  DMSP  SSIES  measures  electron 
temperatures  at  a  relatively  constant  altitude  (840  -  865  km)  as  time  and  position  vary. 
Meanwhile  the  stationary  zenith  antenna  ISR  measures  electron  temperatures  at  a  fixed 
location  above  the  site  as  altitude  and  time  vary.  Major  advantages  of  DMSP  data  are  its 
continuous  measurement  collection  capability  and  better  data  temporal  resolution.  With 
at  least  two  spacecraft  in  operational  at  any  given  time,  data  can  be  collected  for  a 
minimum  of  four  local  times  per  101 -minute  orbit.  DMSP  processes  Te  measurements  in 
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four-second  intervals.  On  the  other  hand,  incoherent  scatter  radars  are  only  typically 
turned  on  for  given  blocks  of  time  to  support  funded  experiments.  ISR  data  profiles  are 
normally  averaged  over  2  to  10  minute  intervals. 

DMSP  and  ISR  each  possess  possible  sources  of  error  imique  to  each  electron 
temperature  measurement  technique.  Rich’s  SSIES  Users  Guide  [1994]  cites  the 
inherent  uncertainty  of  the  SSEES  electron  probe’s  sensitivity  at  ±  50  K.  Since  the  DMSP 
Te  measurement  technique  depends  upon  net  current  collected  by  the  electron  probe, 
unwanted  charge  carriers  that  are  collected,  such  as  photo  and  auroral  electrons,  could 
contaminate  measurements  [Rich,  1994],  Photoelectrons  are  generated  by  any  surface 
exposed  to  sunlight,  while  auroral  electrons  enter  Earth’s  ionosphere  especially  at  times 
of  heightened  geomagnetic  activity.  Rich’s  SSIES  Users  Guide  [1994]  also  asserts  that 
the  contaminating  influence  of  photo  and  auroral  electrons  on  the  total  current  collected 
by  the  EP  will  be  minimal  as  long  as  the  ambient  electron  flux  exceeds  the  photoelectron 
and  auroral  fluxes.  However,  at  lower  ambient  electron  densities,  as  during  solar 
minimum,  an  intense  photoelectron  or  auroral  electron  flux  can  elevate  the  resulting 
electron  temperature  measurements.  As  another  possible  DMSP  Te  source  of  error,  if  the 
spacecraft’s  sensor  potential  (SENPOT)  cannot  be  maintained,  the  SSIES  data  can 
become  unreliable  [Hairston,  private  correspondence,  21  Dec  00] .  For  SENPOT  to 
function  properly,  the  plasma  current  to  the  SENPOT  reference  surface  must  exceed  the 
leakage  and  photoelectron  current  [Rich,  1994]. 

Turning  to  potential  ISR  Te  sources  of  error,  the  number  of  reliable  ISR  Te 
measurements  noticeably  fall  off  around  700  -  800  km  due  to  lower  signal  strength 
arising  from  lower  electron  densities  at  and  above  these  altitudes  [Erickson,  private 
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correspondence,  22  Sep  00].  Consequently,  measurements  especially  at  higher  altitudes 
have  lower  signal-to-noise  ratios  corresponding  to  larger  associated  uncertainties.  In 
addition,  some  of  the  underlying  assumptions  behind  the  ISR  Te  measurement  technique 
may  not  be  met,  perhaps  leading  to  inaccurate  Te  values.  For  example,  plasma  may  not 
demonstrate  a  Maxwellian  velocity  distribution  in  the  presence  of  large  electric  fields,  or 
at  latitudes  and  altitudes  with  low  electron  densities  and  fewer  particle  collisions.  Also, 
Millstone  ISR  data  was  fit  based  on  a  two-ion  (O'*’  and  H'*')  composition  in  the  topside 
ionosphere  [Erickson,  private  correspondence,  14  Sep  00].  Any  deviation  to  the 
assumed  ratio  could  alter  the  ISR  Tg  values.  Furthermore,  recall  that  at  lower 

electron  densities,  the  Debye  length  can  approach  the  radar  wavelength,  generating  an 
ISR  return  shape  exhibiting  a  more  Gaussian  than  dual-humped  shape,  possibly  leading 
to  incorrect  Tg  values. 

POLITE  Campaigns 

The  Plasmaspheric  Observations  of  Light  Ions  in  the  Topside  and  Exosphere 
(POLITE)  observation  campaign  is  a  collaborated  effort  by  ISR  sites  worldwide  to 
maximize  simultaneous  coverage  of  topside  ionospheric  measurements.  Each  of  the  ten 
POLITE  experiments  lasted  from  one  to  three  days  between  winter  1996  and  summer 
2000,  spanning  a  time  window  from  solar  minimum  to  solar  maximum.  Table  3  below 
shows  the  exact  POLITE  experiment  dates  and  times.  Although  the  focus  of  the  POLITE 
initiative  was  to  observe  oxygen,  hydrogen,  and  helium  ions  in  the  topside  ionosphere, 
electron  temperature  and  density  were  measured  as  well.  Since  incoherent  scatter  radars 
for  the  POLITE  campaigns  were  set  to  a  mode  optimizing  returns  from  the  topside 
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ionosphere,  this  provided  an  excellent  opportunity  to  obtain  reliable  topside  electron 


temperatures. 


Tables.  POLITE  1-10 Dates 


POLITE 

Start  Date 

Julian  Day 

Time  (Z) 

Julian  Day 

Time  (Z) 

1 

13-Feb-96 

044 

1553 

045 

1558 

2 

1606 

14-NOV-96 

319 

1639 

3 

3-Jun-97 

154 

1204 

6-Jun-97 

157 

4 

2-Dec-97 

336 

1527 

4-Dec-97 

338 

5 

26-May-98 

146 

1605 

28-May-98 

148 

6 

22-NOV-98 

326 

1422 

24-NOV-98 

7 

8-Oct-99 

281 

1550 

9-Oct-99 

1556 

8 

9-Dec-99 

343 

1610 

lO-Dec-99 

344 

1610 

9 

6- Jan-00 

006 

1304 

7-Jan-OO 

007 

1637 

10 

1-Jul-OO 

183 

1255 

3-Jul-OO 

185 

1956 

Previous  Related  DMSP  Validation  Studies 

Although  numerous  studies  have  been  conducted  validating  ISR  topside 
measurements  with  in  situ,  or  direct,  rocket  measurements,  the  concept  of  comparing 
DMSP  SSIES  measurements  with  ISR  parameters  is  still  in  its  infancy.  A  past  study 
comparing  24  Millstone  Hill  ISR  topside  ion  density  (nj)  measurements  with  DMSP  F8 
and  F9  nj  measurements  was  mentioned  in  Fred  Rich’s  SSIES  User  Guide  [Rich,  1994] . 
Results  of  this  study  revealed  that  in  80  percent  of  the  cases,  DMSP  ni  measurements 
agreed  with  the  ISR  nj  measurement  to  within  20  percent,  which  was  the  limit  of  the  ISR 
accuracy  at  840  km  altitude.  Also,  in  no  case  was  there  a  difference  greater  than  35 
percent.  Rich  and  Sultan  of  Air  Force  Research  Lab’s  (AFRL)  Battlespace  Environment 
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Division  later  completed  a  documented,  more  rigorous  Millstone  Hill  ISR  -  DMSP  ni 
study  by  examining  31  cases  of  near  simultaneous  n;  measurements  taken  between  1989 
and  1991  during  a  solar  cycle  maximum.  The  primary  result  of  that  study  is  that  SSBES 
ion  densities  agree  to  within  a  nine  percent  accuracy  of  the  Millstone  Hill  ISR  ion 
densities  [Sultan  and  Rich,  Draft  Document,  2000] .  A  closing  recommendation  of  that 
research  was  to  pursue  a  DMSP  -  ISR  Te  comparison  study,  hence  driving  this  effort.  To 
date,  no  previous  DMSP  -  ISR  Te  comparison  and  validation  studies  have  been  carried 
out. 
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III.  Methodoloev 


Chapter  Overview 

This  chapter  explains  how  the  DMSP  and  ISR  electron  temperature  (Te)  data  are 
processed,  displayed  and  compared  into  meaningful,  identifiable  results.  First,  a 
description  of  the  raw  DMSP  and  ISR  data  is  given.  The  DMSP  SSIES  provides 
measurements  in  four-second  intervals  without  associated  error,  while  Millstone  Hill’s 
zenith-directed  fixed  antenna  ISR  data  arrives  in  roughly  2-10  minute  intervals  per 
profile  at  three  different  pulse  lengths,  each  with  associated  error.  Meanwhile, 
Sondrestrom’s  ISR  radar  measurements  are  taken  at  only  one  pulse  length  at  13. 5 -minute 
intervals  per  profile.  From  these  data,  DMSP  flyover  times  during  the  POLITE  campaign 
dates  are  identified  over  Millstone  Hill  and  Sondrestrom  ISR  sites.  Then,  using  a 
previous  DMSP  -  ISR  electron  density  paper  as  a  guide  [Rich  and  Sultan,  Draft 
Document,  2000],  we  establish  suitable  criteria  for  designating  near-simultaneous  DMSP 
-  ISR  measurements  for  comparison.  Next  we  employ  averaging,  interpolating, 
extrapolating,  filtering,  and  fit  techniques  to  determine  the  electron  temperature  values 
for  comparison.  Advantages  and  disadvantages  are  then  depicted  for  the  approach  used 
on  each  type  of  data.  Subsequently,  DMSP  data  for  each  conjunction  are  explored  for 
occurrences  of  mid-latitude  ionospheric  features.  In  addition,  cases  comparing  nearly 
concurrent  DMSP  flyover  times  by  different  satellites  are  then  examined.  Both  of  these 
steps  are  designed  to  help  understand  and  justify  DMSP  electron  temperature  behavior. 
The  chapter  concludes  with  a  description  of  how  DMSP  random  noise  over  Millstone 
Hill  is  assessed  and  presented. 
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DMSP  Data  Description 

Dr.  Marc  Hairston  from  the  University  of  Texas  at  Dallas  (UTD)  supplied  the 
DMSP  SSIES  data  corresponding  to  the  POLITE  time  windows  in  chronological  files 
segregated  by  spacecraft  and  orbit.  As  mentioned  earlier,  the  range  of  all  10  POLITE 
campaigns  extended  from  Winter  1996  through  Summer  2000,  thus  covering  solar 
minimum  through  solar  maximum.  SSIES  on  DMSP  spacecraft  F12  -  F15  recorded  the 
measurements  for  this  time  period  (See  Figure  14).  Note  that  F-10  and  F-1 1  data  were 
incomplete  or  not  available  for  this  time  period. 


Figure  14.  Years  DMSP  Spacecraft  Mounted  with  SSIES  are  in  Orbit.  The 
corresponding  stage  of  the  solar  cycle  is  represented  by  the  FI  0.7  cm  Radio  Flux  trend 
[Rich  DMSP  Web  Site.  20001. 


Table  4  summarizes  the  list  of  SSEES  parameters  that  UTD  provided.  The  data 


initially  contained  velocity  components  (except  for  POLITE  10),  and  later  included 


Table  4.  DMSP  Data  Description 


Parameter 

Description 

Time 

Avg.  Measurement  Time  (s) 

Vx 

X  component  of  Ion  Drift  Velocity  (m/s) 

Vy 

Y  component  of  Ion  Drift  Velocity  (m/s) 

Vz 

Z  component  of  Ion  Drift  Velocity  (m/s) 

Density 

Total  Ion  Density  (cm'^) 

MLT 

Magnetic  Local  Time 

MEAT 

Magnetic  Latitude 

glong 

Geographic  Longitude 

glat 

Geographic  Latitude 

Alt 

Spacecraft  Altitude  (km) 

ffach 

Fraction 

frache 

Fraction  He"^ 

ffaco 

Fraction  0"^ 

Ti 

Ion  Temperature  (K) 

Te 

Electron  Temperature  (K) 

SZA 

Solar  Zenith  Angle 

SENPOT 

Sensor  Potential  (V) 

Sensor  Potential  (SENPOT)  and  Solar  Zenith  Angle  (SZA)  to  help  determine  when  the 
satellite  was  sunlit.  Each  complete  orbit  began  and  ended  as  a  spacecraft  crossed  the 
equator  in  its  ascending  node  [Hairston,  private  communication,  24  Aug  00] .  The  data 
were  averaged  into  four  second  bins,  thus  with  each  DMSP  orbit  lasting  around  101 
minutes,  each  file  contained  about  1500  data  points.  Figure  15  shows  a  representative 
DMSP  Te  plot  for  an  entire  orbit. 
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DMSP  Orbit  -  POLITE  4  -  F13 
3  Dec  97  1025-1207Z 
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Figure  15.  Sample  DMSP  Te  Data  for  an  Entire  Orbit 


Before  commencing  the  DMSP  Te  analysis,  AFRL  and  UTD  advised  of  potential 
instrument  related  factors  that  could  distort  some  of  the  Te  data.  These  included  the 
following: 

1.  During  late  1999  and  2000  (POLITEs  7  -  10),  the  F14  sensor  potential  could 
not  remain  at  the  plasma  potential  during  the  simlit  portion  of  the  orbit. 
Consequently,  F14  SSIES  data  is  misrepresented  for  these  daylight  orbit 
sectors  [Hairston,  private  correspondence,  21  Dec  00], 

2.  Due  to  efforts  focused  on  observing  the  Leonid  meteor  shower,  all  POLITE  6 
(22  -  25  Nov  98)  Te  measurements  were  taken  in  30-second  intervals  instead 
of  four  seconds.  [Hairston,  private  correspondence,  21  Dec  00].  This  led  to 
further  smoothing  of  the  Te  data. 

3.  F13  instrument  settings  could  cause  evening  pass  Te  values  to  be 
systematically  high  [Sultan,  private  communication,  26  Sep  00], 
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ISR  Data  Description 

Millstone  Hill  and  Sondestrom  observatories  provided  the  ISR  data  for  the 
POLITE  I  -  10  campaigns.  Dr.  Phil  Erickson  from  Millstone  Hill  supplied  their  ISR 
measurements  in  4  different  pulse  lengths  (410, 640, 1000,  and  2000  ps)  from  both  the 
fixed  and  steerable  antennas.  To  avoid  the  need  to  correct  for  an  ISR  /  DMSP  conjunction 
displacement,  only  data  from  the  fixed  zenith-directed  antenna  (88°  elevation  angle)  were 
used  for  this  study.  Since  most  of  the  640  ps  data  was  missing,  the  three  remaining  ISR 
pulse  lengths  were  used.  Averaged  over  a  smaller  altitude  extent,  the  shorter  ISR  pulse 
lengths  generally  have  smaller  uncertainties  and  offer  finer  altitude  resolution  than  the 
longer  pulse  lengths  [Erickson,  private  correspondence,  22  Sep  00] .  The  410  ps  pulse 
length  provided  the  finest  altitude  resolution  (30  km)  but  only  extended  up  to  75 1  km, 
about  100  km  shy  of  DMSP  altitude.  Meanwhile  the  longer  1000  and  2000  ps  pulse 
lengths  provided  measurements  well  above  DMSP  altitude,  with  vertical  resolutions  of  75 
and  150  km,  respectively.  Table  5  displays  details  of  the  three  Millstone  Hill  ISR  pulse 
lengths.  Note  that  each  altitude  represents  the  center  value  of  the  distance  over  which  the 
pulse  is  averaged.  For  example,  a  1000  ps  Te  measurement  at  874.09  km  with  a  75  km 
vertical  resolution  represents  an  average  Te  value  across  a  150  km  altitude  range  (±75 
km)  centered  upon  874.09  km.  The  ISR  data  for  each  pulse  length  was  integrated  over  a 
one-minute  interval  for  each  profile.  Parameters  used  to  analyze  Millstone  ISR  Te  data  to 
compare  to  SSIES  Te  included  the  following: 

1 .  Time  past  0000  UT  on  the  first  day  of  the  POLITE  campaign 

2.  Altitude  (km) 

3.  Pulse  length 
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4.  Signal-to-Noise  ratio 

5.  Electron  temperature 

6.  Electron  temperature  error 


Table  5.  Millstone  Hill  ISR  Pulse  Length  Details 


Pulse  Length 

41  Ops 

lOOOps 

lOOOps 

Low  Altitude  (km) 

WnSKM 

244.9 

High  Altitude  (km) 

■BHi 

■MQII 

Vertical  Resolution  (km) 

75 

Points  per  Profile 

21 

9 

7 

Altitude  Closest  to  DMSP  (km) 

*751.3 

874.1 

844.1 

*  Highest  Altitude  Available 

As  the  signal-to-noise  (SN)  ratio  decreased  with  altitude,  Dr.  Erickson  cautioned  us  to  be 
wary  of  measurements  with  SN  ratios  <0.1  [Erickson,  private  correspondence,  13  Oct 
2000]. 

Meanwhile,  Mary  McCready  of  SRI  provided  the  Sondrestrom  POLITE 
Campaign  ISR  measurements  at  a  450  ps  pulse  length.  Based  on  McCready’s  description 
of  the  data,  each  Te  altitude  profile  was  generated  from  a  ten-minute  period  where  the 
450  ps  pulse  sampled  the  ionosphere  every  1 1.4  km,  with  the  ISR  return  integrated  over  a 
68  km  altitude  range  parallel  to  the  magnetic  (B)  field.  There  was  a  3.5-minute  time  lag 
between  each  ten-minute  measurement  collection  cycle.  The  altitudes  covered  by  these 
measurements  typically  ranged  from  100  -  1050  km,  although  most  measurements  above 
500  km  were  either  missing  or  contained  error  greater  than  15  percent.  Parameters  used 
to  analyze  Sondrestrom  ISR  Tg  data  to  compare  to  SSIES  Tg  included  the  following; 
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1 .  Measurement  start  time 


2.  Measurement  stop  time 

3.  Altitude  (km) 

4.  Electron  temperature 

5.  Electron  temperature  error 

Extracting  Coinciding  DMSP  -  ISR  Data  Sets 

Incoherent  scatter  radar  (ISR)  is  the  only  ground-based  measurement  technique 
that  can  provide  credible  measurements  of  the  topside  ionospheric  electron  temperature. 
Consequently,  DMSP  flyovers  in  the  vicinity  of  ground-based  radar  sites  offer  the 
opportunity  to  validate  DMSP  SSIES  Te  measurements  with  near-simultaneous  ISR  Te 
values  [Sultan  and  Rich,  Draft  Document,  2000] . 

Overlapping  cases  were  determined  by  filtering  through  the  DMSP  and  Millstone 
Hill  Te  data  sets  to  isolate  cases  where  the  spacecraft  and  ISR  had  measurements  similar 
in  time  and  location  given  certain  cutoff  thresholds.  Based  on  Sultan  and  Rich’s  [2000] 
recent  DMSP  -  ISR  electron  density  paper,  we  filtered  through  the  DMSP  data  set  for 
cases  where  a  DMSP  satellite  flew  within  a  5  degree  latitude-longitude  circle  above  the 
Millstone  Hill  radar  site.  Using  trigonometry,  and  given  Millstone  Hill’s  geographic 
latitude  of  42.62°  N,  this  circle’s  maximum  north  to  south  distance  is  nearly  1112  km, 
and  its  maximum  east  to  west  diameter  is  about  819  km.  A  FORTRAN  program 
searched  for  the  DMSP  measurements  lying  within  such  a  circle.  This  routine  generated 
85  potential  DMSP  overpasses  of  Millstone  Hill.  However,  a  deficiency  of  coexisting 
DMSP  and  ISR  data  left  only  37  such  overpasses  feasible  for  Te  comparison.  Since  the 
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spacecraft  typically  did  not  fly  across  the  full  diameter  of  the  overpass  circle,  of  the  37 
total  cases,  we  foimd  anjwhere  from  six  to  41  points  lying  within  this  circle.  Given  the 
four-second  interval  between  data  points,  this  corresponds  to  the  spacecraft  being  inside 
this  circle  anywhere  from  24  seconds  to  2  minutes  44  seeonds  per  conjunction.  Table  6 
gives  a  summary  of  the  37  cases.  Note  the  case  number  preceding  the  dash  eorresponds  to 
the  POLITE  campaign,  while  the  second  number  denotes  the  chronological  occurrence  of 
the  eonjimction  within  that  campaign.  For  example,  M3-2  denotes  the  second  DMSP 
flyover  within  the  five-degree  circle  above  Millstone  Hill  during  POLITE  Three.  The 
average  DMSP  flyover  time  of  the  points  within  this  five-degree  circle  was  then  used  to 
find  the  eorresponding  concurrent  Millstone  Hill  ISR  data  set.  Shadowing  Rich  and 
Sultan’s  [2000]  cutoff  criteria  established  in  their  DMSP  -  ISR  ion  density  eomparison 
paper,  ISR  Te  measurements  within  ±  30  minutes  of  the  average  DMSP  flyover  time  were 
considered  near  simultaneous  to  DMSP  measurements  within  the  five  degree  cirele.  In 
addition,  at  least  two  sequential  altitude  profiles  within  the  ±  30-minute  time  window 
were  required  to  declare  the  ISR  case  as  a  candidate  for  comparison  with  DMSP. 
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Table  6.  Case  Summary  for  DMSP  -  Millstone  ISR  Te  Comparisons 


Case# 

POLITE 

Satellite 

Points  Inside  5®  Circle 

Ml-l 

1 

96045 

0144 

F12 

41 

Ml-2 

1 

96045 

1047 

F13 

29 

M2-1 

2 

96317 

1501 

F12 

31 

M2-2 

2 

96318 

1449 

F12 

41 

M2-3 

2 

96319 

1436 

F12 

36 

M3-1 

3 

97154 

1427 

i  F12 

18 

M3-2 

3 

97155 

0154 

!  F12 

13 

M3-3 

3 

97155 

mm 

29 

M3-4 

3 

97155 

■39 

22 

M3-5 

3 

97156 

0142 

F12 

39 

M4-1 

4 

97337 

1103 

F13 

40 

M4-2 

4 

97337 

1400 

F14 

39 

M4-3 

4 

97337 

1423 

F12 

26 

M4-4 

4 

97338 

1051 

F13 

26 

M4-5 

4 

97338 

1347 

F14 

38 

M5-1 

5 

0055 

F14 

39 

M5-2 

5 

1053 

F13 

27 

M5-3 

5 

98147 

2221 

F13 

23 

M5-4 

5 

98148 

1042 

F14 

34 

M6-1 

6 

98327 

0038 

F14 

11 

M6-2 

6 

98327 

1114 

■391 

40 

M6-3 

6 

98327 

1414 

BHHB 

37 

M7-1 

7 

99282 

41 

M7-2 

7 

99282 

0114 

27 

M7-3 

7 

99282 

1100 

F13 

24 

M8-1 

8 

99344 

0121 

F14 

15 

M8-2 

8 

99344 

1134 

F13 

24 

M8-3 

8 

99344 

22 

M9-1 

9 

00006 

■391 

28 

MlO-1 

10 

00183 

1348 

F14 

29 

MlO-2 

10 

00183 

1434 

F15 

39 

MlO-3 

10 

00183 

2208 

F13 

40 

MlO-4 

10 

00184 

0115 

F14 

22 

MlO-5 

10 

00184 

1421 

F15 

12 

Ml  0-6 

10 

00184 

2156 

F13 

21 

Ml  0-7 

10 

00185 

0102 

F14 

40 

Ml  0-8 

10 

00185 

0146 

F15 

32 

The  same  five-degree  circle  search  routine  was  used  to  identify  DMSP  overpasses 
above  the  Sondrestrom  radar  site.  This  yielded  169  potential  cases,  of  which  56  had 
corresponding  ISR  data  available.  To  account  for  the  short  scale  length  and  large 
variations  of  ionospheric  phenomena  typical  at  auroral  latitudes,  the  five-degree  circle 
was  reduced  to  one  degree.  Since  Sondrestrom’s  ISR  signal  points  up  the  magnetic  (B) 
field  line  with  an  approximate  elevation  angle  of  80°,  this  one-degree  circle  was  slightly 
offset  fi'om  being  directly  above  Sondrestrom’s  location  to  account  for  the  departure  of 
the  beam  fi'om  vertical.  The  smaller  circle  resulted  in  15  DMSP  flyovers,  seven  of  which 
had  corresponding  ISR  data.  Table  7  shows  a  breakdown  of  these  seven  DMSP  - 
Sondrestrom  conjunctions.  Typically  three  to  six  data  points  were  found  within  the  circle 


Table  7.  Case  Summary  for  DMSP  -  Sondrestrom  ISR  Te  Comparisons 


Case# 

POLITE 

Spacecraft 

YearDay 

Time  (Z) 

Pts.  Inside  1®  Circle 

S2-1 

2 

F12 

96316 

2324 

6 

S3-1 

3 

F14 

97156 

2222 

6 

3 

97157 

1018 

5 

S5-1 

5 

98146 

1320 

5 

S64 

6 

98328 

2238 

6 

S8-1 

8 

F13 

99343 

1957 

3 

S9-1 

9 

F13 

00007 

1039 

4 

for  a  given  flyover.  The  variability  of  the  DMSP  Te  data  even  within  this  one-degree 
circle  can  still  be  rather  large  to  analyze  confidently.  Figure  16  shows  such  an  example 
of  a  Sondrestrom  case. 


46 


F13  DMSP  Te  ±  3  min  of  Sondrestrom 


Geographic  Latitude  Avg.  Te  within  1“:  2100K 
CT  =  504.65  K 


Figure  16.  DMSP  (FI 3)  Te  Data  Near  a  Sondrestrom  Flyover  on  7  Jan  00. 
Notice  the  large  Te  variability  over  a  short  distance  typical  at  auroral  latitudes 
(About  a  1000  K  fluctuation  within  1°  of  Sondrestrom  (Bold  Line)). 


In  Figure  16,  the  standard  deviation  (c)  of  504.65  K  within  the  one-degree  circle  of 
Sondrestrom  more  than  doubles  the  average  cr  of  242.93  K  calculated  for  the  37  cases  of 
DMSP  Te  measurements  within  the  5  degree  circle  of  Millstone  Hill.  Consequently  it  is 
much  more  difficult  to  distinguish  between  noise  related  and  physics  related  Te  behavior 
at  high  latitudes  than  at  mid-latitudes. 


Calculating  DMSP  Te  Comparison  Values 

Once  the  37  DMSP  -  Millstone  ISR  conjunctions  were  identified,  average  Te 
values  were  calculated  to  determine  the  final  Te  comparison  values.  The  DMSP  values 
used  for  comparison  merely  involved  averaging  the  Te  measurements  that  fell  within  the 
five-degree  circle  for  each  case.  Table  8  below  shows  mean  DMSP  electron  temperatures 
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for  each  case.  Note  that  the  Te  values  obtained  for  Cases  M8-3,  MlO-1,  MlO-4,  and 
Ml 0-7  are  unrealistic,  most  likely  due  to  a  probable  SENPOT  error  to  be  discussed  in 
Chapter  IV.  Case-by-case  plots  of  DMSP  Te  measurements  over  Millstone  Hill  are 
available  as  a  supplementary  document  upon  request. 


Table  8.  Mean  DMSP  Te  Value  Case  Summary 


Case# 

Satellite 

Mean  Te  (K) 

Ml-1 

F14 

2974.5 

Ml-2 

3455.0 

M2-1 

mm 

4988.4 

3693.5 

M2-2 

F12 

4522.9 

2486.8 

M2-3 

4716.7 

■B 

2580.7 

M3-1 

F12 

4356.7 

M7-3 

F13 

3505.8 

M3-2 

F12 

3211.5 

M8-1 

mm 

3129.3 

M3-3 

F13 

4071.7 

M8-2 

■B 

3005.8 

M3-4 

3856.4 

’M8-3 

F12 

5592.7 

M3-5 

3412.8 

M94 

F13 

3318.6 

M4-1 

mm 

2960.5 

‘MlO-l 

■B 

811.0 

M4-2 

4146.2 

MlO-2 

■B 

4144.6 

M4-3 

F12 

4223.8 

M10<3 

■B 

3469.0 

M4-4 

F13 

2903.1 

‘MlO-4 

bh 

779.1 

M4-5 

F14 

4404.2 

Ml  0-5 

F15 

3895.0 

M5-1 

F14 

2702.6 

MlO-6 

F13 

3305.7 

M5-2 

F13 

4084.4 

’MlO-7 

F14 

772.0 

M5-3 

F13 

3821.7 

MlO-8 

F15 

2839.4 

M5-4 

F14 

2822.4 

Notes:  '  Probable  SENPOT  error. 


Similarly,  the  mean  DMSP  electron  temperature  for  the  seven  possible  Sondrestrom  cases 
was  calculated  by  averaging  the  data  points  within  the  one-degree  circle  (“up  B”)  the 
same  way;  Table  9  reveals  these  values. 
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Table  9.  Sondrestrom  Mean  DMSP  Te  Value  Case  Summary 


Case  # 

Satellite 

Mean  T*  (K) 

S2-1 

F12 

S3-1 

F14 

2726.7 

S3-2 

F13 

S5-1 

F14 

3824.0 

S6-1 

F14 

4320.0 

S8-1 

F13 

4036.7 

S9-1 

F13 

2160.0 

Calculating  Millstone  ISR  Comparison  Values 

Determining  the  average  ISR  Te  values  at  DMSP  altitude  for  each  case  involved  a 
combination  of  averaging,  fitting,  interpolation,  and  extrapolation  techniques  depending 
on  the  ISR  mode.  On  average,  about  15  profiles  from  the  1000  ps  and  2000  ps  data  sets 
fell  within  the  ±  30  minutes  of  DMSP  flyover  time,  compared  to  typically  three  or  four 
profiles  for  each  410  ps  data  set.  The  maximum  altitude  and  vertical  resolution  of  each 
pulse  length’s  data  helped  dictate  the  method  used  to  compute  an  average  Te  value  at 
DMSP  altitude  (refer  back  to  Table  5). 

Since  the  410  ps  profiles  do  not  extend  up  to  DMSP  orbit  altitude  (840  -  865  km), 
extrapolation  using  a  simple  curve  fitting  routine  was  required  to  estimate  Te  values. 

Note  that  associated  Te  error  was  not  extrapolated  with  the  Te  measurement  due  to  the 
complexity  of  such  a  procedure.  Due  to  a  small  signal-to-noise  (SN)  ratio  and 
corresponding  larger  uncertainty  at  higher  altitudes,  a  data  “filter”  was  applied  to  boost 
prospects  for  a  more  meaningful  comparison  to  DMSP  Te  values.  A  FORTRAN  program 
searched  through  the  410  ps  profiles  already  within  ±  30  minutes  of  DMSP  flyover  time 
for  Te  measurements  that  met  the  following  criteria: 
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1.  SN  ratio  >0.1 


2.  Te  error  <10  percent  relative  error 

3.  Te  +  Te  Error  Te  at  the  next  lowest  altitude  (i.e.,  Te  should  not  decrease  with 
altitude.) 

In  addition,  to  carry  on  with  the  extrapolation,  there  had  to  be  at  least  two  resulting 
“filtered”  Te  vs.  altitude  profiles,  each  with  at  least  one  data  point  above  450  km.  Figures 
17  and  18  represent  examples  of  “unfiltered”  and  “filtered”  410  ps  profiles,  respectively, 
for  a  given  conjunction.  Filtering  eliminated  most  of  the  higher  altitude  points  with  large 
error,  leaving  behind  more  reliable  profiles  for  extrapolation. 


Figure  18.  Filtered  410  |j,s  Te  Data 


A  simple  curve  fit  using  Table  Curve2D  software  provided  an  extrapolated  Te  value  to 
DMSP  altitude  for  each  “filtered”  Te  profile.  For  several  cases,  some  of  the  lower 
altitude  points,  especially  below  200  km,  were  removed  to  provide  the  most 
representative  fit  through  the  majority  of  the  topside  data.  If  all  the  Te  profiles  for  a 
given  conjunction  were  close  enough  together  upon  visual  inspection,  only  one  curve  fit 
and  extrapolation  was  performed  to  determine  the  Te  comparison  value.  Otherwise,  a 
curve  fit  was  accomplished  on  each  profile,  with  the  resulting  extrapolated  values 
averaged  to  produce  the  Te  comparison  value.  Thirty-five  of  the  37  conjunctions  met  all 
of  the  criteria  to  compute  extrapolated  Te  comparison  values  from  the  410  ps  profiles. 
Figure  19  shows  an  example  of  the  curve  fitting  and  extrapolation  process  on  one 
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“filtered”  410  fxs  profile.  Appendix  B  contains  a  list  of  averaged  410  ps  extrapolated  Te 
comparison  values,  while  the  plots  that  generated  these  values  are  available  upon  request. 


Figure  19.  Extrapolation  Using  a  Simple  Curve  Fit  of  One  of  Two  “Filtered”  410 
ps  Profiles  Available  for  a  Conjunction.  The  profile  with  the  error  bars  is  the  one 
being  fit.  Error  is  not  extrapolated  with  the  Te  measurement.  Notice  the  data  point 
at  1 82  km  is  not  included  to  aid  in  smoothing  the  fit  through  the  rest  of  the  data 
points.  Extrapolated  Te  value  at  DMSP  altitude  (854  km)  is  2475.1  IK 


Meanwhile,  the  1000  and  2000  ps  ISR  Te  profiles  contained  the  DMSP  altitude 
within  its  altitude  limits,  thus  averaging  pre-existing  measurements  close  to  DMSP 
altitude  to  generate  the  Te  value  for  DMSP  comparison  was  possible.  Just  as  with  the  410 
ps  data,  a  representative  data  filter  was  applied  to  1000  and  2000  ps  ISR  data  before  any 
comparison  Te  values  were  calculated.  However,  this  time  the  Te  data  filter  was  adjusted 
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to  handle  the  longer  pulse  length  and  higher  altitude  measurements  by  searching  for  Tg 
data  meeting  the  following  criteria: 


1 .  Tg  error  <15  percent  relative  error 

2.  Tg  +  Tg  Error  Tg  at  the  next  lowest  altitude  (same  as  with  41  Ops  filter) 

3.  SN  ratio  threshold  ratio  not  used 

These  filter  thresholds  were  relaxed  slightly  compared  to  the  410  ps  Tg  cutoff  criteria 
mainly  to  accommodate  reasonable  Tg  measurements  close  to  DMSP  altitude  with  errors 
between  10-15  percent  and/or  a  SN  ratio  just  below  0.1.  Such  values  would  have  been 
dropped  imder  the  410  ps  filter  thresholds.  To  support  this  assertion.  Table  10  shows  a 
1000  ps  profile  collected  during  POLITE  6.  In  this  case,  the  data  points  at  949  and  1023 
km  (bolded  and  shaded)  would  be  removed  by  both  filters,  while  the  data  points  at  799 
and  874  km  (bolded)  would  survive  the  adjusted  filter  threshold,  but  would  be  removed 


Table  10.  Filter  Application  on  POLITE  6  Millstone  Hill  1000  ps  ISR  Profile 


(14.44Z  23  Nov  98) 


SN  Ratio 

Te  Error 

Rel.  %  Error 

14.44 

424.67 

1.899 

2342.2 

17.8 

0.76% 

14.44 

499.57 

0.767 

2603.8 

21.3 

0.82% 

14.44 

574.47 

0.403 

2637.1 

44.4 

1.68% 

IS9II 

649.38 

0.223 

118.9 

4.21% 

724.28 

0.135 

198.3 

6.80% 

799.18 

0.099 

IH9i 

301.6 

9.98% 

14.44 

874.09 

3158.2 

330.0 

10.45% 

14.44 

14.44 

i 
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by  the  original  410  |^s  filter  for  failing  to  meet  the  <10  percent  error  and  >0.1  SN  ratio 
requirements.  Figure  20  contains  a  plot  for  this  same  profile  to  show  that  the  874  and 
799  km  data  points  that  siirvive  the  less  restrictive  filter  appear  to  be  part  of  the  profile 
trend,  thus  it  seems  reasonable  to  retain  such  data.  Had  these  points  been  subject  to  the 
tighter  filter,  averaging  and/or  interpolating  this  profile  would  not  have  been  possible. 


14.44Z  23  Nov  98  Millstone  Hill  POLITE  6 


1100-1 

Filter  Analysis  on  a  1000  ps  Profile 
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Figure  20.  Impact  of  the  Two  Filters  on  a  Millstone  1000  ps  POLITE  6  Profile.  Notice 
how  the  more  restrictive  filter  (<  10  percent  error  and  >0.1  SN  ratio  thresholds) 
eliminates  two  points  that  seem  to  be  part  of  the  profile  trend. 


On  average,  ISR  electron  densities  gradually  increase  to  be  about  one  order  of 
magnitude  greater  during  POLITE  10  (solar  maximum)  than  during  POLITE  1  (solar 
minimum)  [Keyser,  private  correspondence,  12  Feb  01] .  Higher  electron  densities  imply 
a  stronger  ISR  signal,  thus  a  higher  SN  ratio  and  lower  absolute  error  for  an  ISR 
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measurement  at  a  given  altitude.  This  Ue  versus  solar  cycle  trend  perhaps  explains  why 
more  ISR  Te  measmements  survive  the  filtering  process  towards  solar  maximum. 

The  averaging  technique  for  both  1000  and  2000  ps  pulse  lengths  involved 
extracting  and  calculating  the  mean  of  the  “filtered”  Te  measurements  closest  to  DMSP 
altitude.  These  measurements  occurred  at  874  km  for  the  1000  ps  data,  and  at  844  km  for 
the  2000  ps  data,  while  actual  DMSP  flyover  altitudes  ranged  from  841  -  863  km.  Since 
ISR  Te  measurements  closest  to  DMSP  altitude  are  already  averaged  over  a  vertical  range 
including  the  DMSP  altitude  in  the  first  place,  differences  associated  with  averaging  at 
844  or  874  km  versus  at  the  exact  DMSP  altitude  should  be  small. 

To  defend  this  assumption,  linear  interpolation  to  DMSP  altitude  and  then 
averaging  the  interpolated  values  was  performed  on  the  1000  ps  data  to  see  just  how 
different  the  Te  comparison  results  were  from  averaging  at  874.09  km.  This  procedure 
was  only  carried  out  on  the  1000  ps  data  because  the  2000  ps  upper  altitude  bounding 
DMSP  altitude  (993  km)  required  for  linear  interpolation  was  usually  missing  or  did  not 
make  the  filter  due  to  relative  Te  error  >15  percent.  Linear  interpolation  was  chosen 
instead  of  a  spline  or  polynomial  fit  because  the  Te  gradient  in  the  75  km  interval 
containing  DMSP  altitude  is  typically  quite  small  [Erickson,  private  correspondence,  13 
Oct  00].  For  a  given  1000  ps  profile,  linear  interpolation  to  DMSP  altitude  was  only 
possible  as  long  as  the  two  data  points  sandwiching  the  DMSP  altitude  (874  km  and  799 
km)  survived  the  filtering  process.  The  average  of  the  resulting  linearly  interpolated 
values,  as  long  as  there  were  at  least  two  “filtered”  profiles  available  per  conjunction, 
produced  the  final  Millstone  Te  value  compared  against  DMSP.  Table  1 1  shows  a 
comparison  of  linearly  interpolated  Te  values  to  DMSP  altitude  with  Te  values  averaged 
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at  874.09  km.  Note  that  of  the  37  total  cases,  the  16  cases  in  the  table  represent  those  that 
survived  the  filtering  process,  and  thus  were  suitable  for  interpolation.  In  all  16  of  these 
cases,  the  averaged  Te  value  closest  to  DMSP  altitude  exceeded  the  linearly  interpolated 
Te  average  at  the  exact  DMSP  flyover  altitude  (874.09  km)  by  an  average  of  52.80  K. 

This  average  difference  is  about  five  times  less  than  the  average  absolute  Te  error 
associated  with  the  Te  values  themselves.  Thus  the  difference  between  the  two 
procedures  is  relatively  insignificant. 


Table  11.  Comparison  of  Averaging  vs.  Linear  Interpolation  Techniques 


Case# 

lOOO^s  Te  Avg 

Avg. 

Error 

*  lOOOns  Te  Lint 

Avg. 

Error 

M5-1 

2139.20 

237.86 

2133.27 

211.09 

5.93 

3028.57 

314.05 

2994.99 

289.61 

2266.20 

269.96 

2230.20 

246.20 

Bsai 

327.80 

3027.40 

103.74 

mam 

308.03 

2420.48 

291.46 

54.52 

M7-2 

2456.94 

300.83 

2366.04 

M8-3 

3014.86 

353.60 

3011.52 

M9-1 

2766.80 

242.70 

2756.64 

210.90 

10.16 

3576.95 

373.36 

3464.35 

112.6 

3637.09 

419.75 

3511.77 

125.32 

3209.01 

186.32 

3151.90 

173.19 

57.11 

2376.51 

156.68 

2330.30 

139.76 

46.21 

3166.29 

3154.48 

254.37 

11.81 

2868.91 

164.75 

2830.24 

148.11 

3S.67 

2437.31 

2383.69 

2170.73 

2109.39 

Averages 

2795.09 

265.43 

2742.29 

235.16 

52.80 

Note:  ^  “Lint”  =  Linear  Interpolated 
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Determining  Millstone  ISR  Uncertainty 

A  representative  Millstone  Hill  ISR  uncertainty  provides  one  of  the  benchmarks 
against  which  the  overall  significance  of  the  DMSP  -  ISR  percent  difference  can  later  be 
compared.  Each  ISR  measurement  we  received  included  a  Te  value  and  an  associated 
absolute  error,  or  absolute  uncertainty  (ATe).  Since  Te  values  are  averaged  closest  to 
DMSP  altitude  to  determine  a  mean  Te  comparison  value  for  the  1000  and  2000  ps 
measurements,  standard  error  propagation  methods  are  employed  to  express  a 
conjunction’s  Te  comparison  uncertainty  as  the  average  of  individual  uncertainties,  or 
explicitly: 

(11) 

n 

where 

(Are)„yg  =  Average  Absolute  Te  Uncertainty 

n  =  Number  of  ATe  Values 
ATe/  =  “i**”’  Absolute  Te  Uncertainty  Value 

Using  this  method,  average  ISR  Te  uncertainties  can  be  calculated  by  conjunction,  pulse 
length,  and  across  all  pulse  lengths.  Note  that  since  the  410  ps  ISR  measurements  were 
extrapolated  to  DMSP  altitude  without  incorporating  error,  the  overall  Millstone  Hill  ISR 
uncertainty  presented  in  Chapter  IV  is  only  based  upon  available  1000  and  2000  ps  Te 
comparison  values. 
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Calculating  Sondrestrom  ISR  Comparison  Values 

Since  most  Sondrestrom  ISR  Te  data  from  over  500  km  were  either  missing  or 
contained  error  greater  than  15  percent,  the  smoother  lower  altitude  data  profiles  were 
extrapolated  to  DMSP  altitude.  Like  with  the  Millstone  ISR  data,  at  least  two 
consecutive  profiles  within  ±30  minutes  of  DMSP  flyover  time  had  to  exist  for  a 
conjunction  to  be  considered  for  extrapolation,  fri  addition,  a  relative  error  threshold  of 
15  percent  was  established  as  a  filtering  criterion.  Due  to  the  variable  structure  and 
limited  number  of  “filtered”  measurements  available,  especially  above  500  km, 
determining  an  extrapolated  average  Te  at  DMSP  altitude  proved  to  be  quite  complicated. 

DMSP  -  ISR  Te  Comparison  Statistics 

Once  DMSP  and  Millstone  Hill  ISR  Te  comparison  values  were  obtained  using 
the  above  methods,  percent  differences  between  each  conjunction’s  DMSP  Te  and  each 
ISR  pulse  length’s  Te  value  were  calculated.  In  each  case,  the  percent  difference  was 
computed  according  to: 


Percent  Difference  = 


(DMSP  Te  -ISR  Te) 

^  X  100 

ISRTe 


(12) 


where  DMSP  Te  and  ISR  Te  represent  DMSP  and  ISR  Te  values,  respectively.  The 
average  and  standard  deviation  of  the  percent  differences  for  all  cases  were  then 
calculated  to  determine  an  overall  DMSP  percent  difference  and  standard  deviation 
compared  to  the  ISR  uncertainty  and  the  quoted  SSIES  Te  measurement  uncertainty  of 
ten  percent  [Rich,  1994].  Average  DMSP  Te  -  ISR  Te  percent  differences  are 
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summarized  by  pulse  length  in  Chapter  IV,  and  presented  by  conjunction  and  pulse 
length  in  Appendix  A. 

Case  Study  Methodology 

Of  the  37  Millstone  Hill  conjunctions,  a  number  of  concurrent  DMSP  Te  and 
electron  density  (ne)  profiles  were  examined  to  better  understand  SSIES  Te  data 
consistency,  SSIES  instrument  performance,  and  to  identify  ionospheric  phenomena  that 
the  data  revealed.  Instances  demonstrating  SSIES  data  consistency  involved  comparing 
Te  data  measured  by  different  spacecraft  flying  through  the  same  five-degree  circle  over 
Millstone  Hill  within  tens  of  minutes  of  each  other.  Occmrences  of  DMSP  Te  profile 
behavior  outside  physical  and  noise  level  expectations  were  inspected  for  possible 
instrument  related  anomalies.  Meanwhile  the  structure  of  simultaneous  Te  and  ne  profiles 
were  compared  to  look  for  such  mid-latitude  features  as  light  ion  troughs,  mid-latitude 
electron  density  troughs,  and  subauroral  ion  drifts  (SAID).  See  Table  18  in  the  “Results” 
Chapter  for  the  case  summary  breakdown. 

Determining  DMSP  Te  Noise  Level 

As  mentioned  in  the  introduction,  one  of  the  goals  of  this  research  was  to  assess 
the  DMSP  Te  random  noise  level,  as  this  is  crucial  information  for  space  weather 
modeling  efforts  such  as  GAIM.  The  noise  level  quantification  procedure  based  on  a 
linear  Te  variation  with  latitude  represents  a  first  attempt  at  this  objective.  DMSP  Te  data 
samples  were  restricted  to  the  same  37  cases  within  the  five-degree  circle  above 
Millstone  Hill  that  were  used  for  the  ISR  comparison.  Consequently,  all  situations 
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analyzed  were  located  at  mid-latitudes  in  the  vicinity  of  42.6  N  (53.2  Geomagnetic 
North).  Note  that  approximating  a  linear  DSMP  Te  behavior  for  the  Sondrestrom 
conjunctions  would  have  been  unreasonable  due  to  the  short  scale  lengths  of  ionospheric 
behavior  at  high  latitudes.  The  final  statistic  obtained  to  gauge  a  data  set’s  noise  level 
was  the  standard  deviation  as  fraction  of  the  filtered  mean  electron  temperature  (SDFM). 
The  steps  used  to  obtain  the  SDFM  included  the  following: 

1 .  Perform  a  linear  regression  on  the  raw  data  of  form:  y  =  a  +  bx. 

2.  Determine  residuals  by  subtracting  the  fit  from  the  raw  data. 

3.  Compute  the  standard  deviation  of  the  residuals  (crro). 

4.  Remove  data  points  with  residuals  outside  of  two  standard  deviations  (2  a  ro) 
of  the  residual  distribution.  The  resulting  data  set  is  termed  the  “filtered”  data 
set. 

5.  Compute  a  mean  Te  ( Te^  )  of  the  filtered  data. 

6.  Perform  a  linear  regression  on  the  filtered  data. 

7.  Determine  residuals  from  the  filtered  data. 

8.  Compute  the  standard  deviation  of  these  residuals  (ar{). 

9.  Divide  crri  by  Te,  to  obtain  the  standard  deviation  as  a  fi'action  of  the  filtered 
mean  Te.  Thus  SDFM  = 

Based  on  the  assumption  that  Te  varies  linearly  with  latitude,  the  simple  linear  fit 
was  applied  with  the  intent  to  remove  physical  trends  in  the  Te  data.  Thus  any  Te 
variations  around  the  fit  line  were  considered  random  noise.  Points  lying  outside  the  2a 
boundary  (Step  4)  were  assumed  to  have  at  least  a  95  percent  probability  of  being  “bad,” 
and  were  removed.  This  procedure  is  a  simplified  version  of  Chauvenet’s  Criterion  for 
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rejection  of  data,  where  a  data  point  is  rejected  if  the  probability  of  its  occurrence  (based 
on  a  Gaussian  distribution)  multiplied  by  the  number  of  data  points  in  the  set  is  50 
percent  or  less /Ma^/ue.s'en,  1997]. 

The  fact  that  Te  does  not  truly  behave  linearly  with  respect  to  latitude  poses  the 
risk  of  poor  line  fits,  which  in  some  cases  could  filter  out  Te  data  points  that  appear 
reasonable.  Figures  21-24  below  graphically  step  through  an  example  of  the  DMSP  Te 
noise  level  determination  process  to  obtain  the  SDFM.  All  SDFM  results  are  presented 
in  Table  18  in  Chapter  IV. 


F13  DMSP  Te  Within  5°  Circle  of  Millstone  Hill 


Geographic  Latitude  Mean  Te:  3469.0  K 

ti  of  Residuals:  109.92  K 


Figure  2 1 .  Initial  Linear  Fit.  One  point  lies  outside  of  2<j  boundary. 
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Figure  22.  Plot  of  Residuals  Based  on  Linear  Fit.  One  point  lies  outside  of 
2a,  thus  is  removed  before  the  new  linear  fit. 


Figure  23.  Linear  Fit  of  Filtered  Data 
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Residuals  (K) 


"Filtered"  F13  DMSP  Te  Residuals  Within  5"  Circle  of  Millstone  Hill 
POLITE  10  -22.1 3Z  (2208Z)  1  Jul  00  Flyover 

200 

100 

0 

-100 

-200 

-300  J 

Geographic  Latitude 

CT  of  Residuals:  102.69  K 
a  as  Fraction  of  Te  Mean:  .030 

Figure  24.  Plot  of  “Filtered”  Residuals 
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IV.  Results  and  Analysis 


Chapter  Overview 

Averaged  DMSP  Te  values  for  all  37  Millstone  flyovers  and  six  Sondrestrom 
flyovers  were  compared  to  corresponding  ISR  filtered  Te  estimates  at  DMSP  altitudes. 
Three  DMSP  FI 4  Millstone  flyovers  during  POLITE  10  and  an  F12  during  POLITE  8 
produced  unrealistic  values  due  to  a  SSIES  SENPOT  malfunction  and  were  not  included 
in  the  result  tally.  Also,  due  to  missing  data,  or  not  enough  data  surviving  the  filtering 
process  described  in  the  previous  chapter,  ISR  Te  estimates  from  all  pulse  lengths  were 
not  obtainable  for  every  conjunction,  but  there  was  at  least  one  ISR  Te  estimate  available 
per  each  conjunction.  When  more  than  one  Millstone  Hill  ISR  pulse  length  Te  estimate 
was  available,  their  values  lay  well  within  each  other’s  relative  error. 

The  behavior  of  some  DMSP  Te  profiles  showed  indications  of  possible 
instrument  and  telemetry  related  anomalies  that  led  to  extended  segments  of  unrealistic  or 
flat-lined  electron  temperatiure  values.  Nonetheless,  the  37  cases  of  near  concurrent 
Millstone  Hill  ISR  and  DMSP  electron  temperature  measurements  for  the  POLITE  1-10 
campaigns  reveal  that  DMSP  Te  values  consistently  exceed  the  complementing  raw, 
filtered,  and  averaged  ISR  measurements.  The  average  percent  difference  between 
DMSP  and  Millstone  ISR  Te  values  was  nearly  25  percent,  exceeding  ISR  uncertainty  by 
nearly  three-to-one.  Corresponding  DMSP  random  uncertainties  found  later  in  this  study 
(about  four  percent)  do  not  substantially  alter  this  conclusion.  The  margin  between 
DMSP  and  ISR  Te  measurements  decreases  with  time  fi:om  solar  minimum  towards  solar 
maximum,  partially  due  to  a  decreased  photoelectron  influence  on  DMSP  Te 
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measurements  towards  solar  maximum.  Although  high  latitude  physical  variations  over 
short  scale  lengths  can  complicate  the  analysis,  a  brief  look  at  the  seven  DMSP  - 
Sondrestrom  ISR  Te  comparisons  show  similar  results.  The  nine  Millstone  instances 
where  an  ISR  Te  average  exceeded  its  DMSP  counterpart  happen  to  occur  when  there  is 
an  argument  for  faulty  SEESS  performance  or  a  large  uncertainty  in  the  ISR  averaging 
process.  Plots  and  tables  supporting  these  results  are  found  in  this  chapter  and  in  the 
appendices.  Meanwhile  corresponding  DMSP  electron  density  measurements  are 
repeatedly  lower  than  the  ISR  average,  verifying  the  expected  inverse  relationship 
between  electron  temperature  and  density  for  a  fixed  amount  of  energy  [Keyser,  Thesis 
Draft.  2001 J. 

Next,  the  variability  of  DMSP  measured  Te  data  over  Millstone  Hill  is  assessed  to 
provide  a  first  guess  of  a  mid-latitude  random  noise  level  to  be  considered  by  GAIM’s 
Kalman  filter.  In  general,  the  average  DMSP  electron  temperature  random  noise  level  for 
the  Millstone  Hill  conjunctions  is  under  four  percent,  which  is  well  within  the  electron 
probe’s  published  accuracy  of  10  percent  [Rich,  1994].  Details  of  each  case’s  random 
noise  calculation  are  presented  in  a  table.  Noise  level  is  not  determined  for  the 
Sondrestrom  conjunctions  in  order  to  avoid  complications  caused  by  the  short  scale 
lengths  of  high  latitude  ionospheric  events. 

As  part  of  the  DMSP  Te  validation  process,  case  studies  are  analyzed  to  support 
specific  electron  temperature  behavior  and  structure.  When  comparing  Te  profiles  from 
different  spacecraft  with  similar  Millstone  Hill  flyover  times,  electron  temperature 
measurements  are  consistent  within  the  aforementioned  random  noise  level.  Finally,  the 
chapter  closes  with  a  look  at  some  cases  showing  indications  of  mid-latitude  ionospheric 
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features  sueh  as  sub-auroral  ion  drifts,  light  ion  troughs,  and  mid  latitude  electron  density 
(lie)  troughs. 

Millstone  Conjunctions 

Based  on  the  filtering,  averaging,  and  extrapolation  techniques  introduced  in  the 
previous  chapter,  at  least  one  Millstone  Hill  ISR  Te  estimate  could  be  calculated  per 
conjunction.  Of  the  37  Millstone  Hill  conjunctions,  thirty-five  410  ps,  twenty  1000  ps, 
and  twenty-six  2000  ps  filtered  ISR  Tg  values  could  be  calculated.  Table  12  shows  a  case 
breakdown  of  these  values  and  associated  error.  Recall  that  the  1000  and  2000  ps  pulse 
lengths  are  averaged  to  874  and  844  km,  respectively  and  that  410  ps  Tg  values  do  not 
include  an  average  error  due  to  the  complexity  of  computing  an  extrapolated  average 
error.  The  average  Millstone  Hill  ISR  uncertainty  based  on  all  available  averaged  1000 
and  2000  ps  values  was  8.42  percent.  Only  filtered  data  were  used  to  determine  this  ISR 
uncertainty  in  order  to  generate  the  most  representative  ISR  Tg  variation  at  DMSP 
altitude. 

Due  in  part  to  a  higher  rig  toward  solar  maximum,  more  ISR  Tg  data  survived  the 
filtering  process  from  the  more  recent  POLITE  campaigns  than  from  the  earlier  ones. 
Also,  when  more  than  one  Millstone  Hill  ISR  pulse  length  Tg  estimate  was  available  for  a 
given  conjunction,  most  of  these  values  lay  well  within  each  other’s  average  relative 
error.  This  demonstrates  a  consistency  amongst  the  three  pulse  length  Tg  estimates  that 
were  compared  against  DMSP  Tg  values.  The  values  that  did  not  fall  within  each  pulse 
length’s  relative  error  were  typically  either  based  on  only  two  profiles  or  had  SN  ratios 
less  than  0.1. 
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Table  12.  Millstone  Hill  ISR  Te  Estimates  at  DMSP  Altitude 


Case# 

Sat. 

DMSP  Alt.  (km) 

410  ps  Te 

1000  ps  Te 

Avg  Error 

2000  ps  Te 

Avg  Error 

Ml-1 

F12 

861.87 

2207.91 

No  Data 

UtfAHW 

154.15 

Ml-2 

I3CT 

855.62 

* 

No  Data 

'3110.90 

222.40 

M24 

846.84 

2957.98 

* 

* 

M2-2 

ijfci 

847.57 

2600.31 

♦ 

* 

M2-3 

Idfel 

846.77 

2890.01 

* 

3473.57 

M3-1 

F12 

860.82 

3088.92 

* 

2905.23 

M3-2 

F12 

862.35 

2905.32 

239.75 

* 

M3-3 

F13 

850.13 

2530.55 

* 

* 

M3-4 

F13 

855.62 

* 

* 

2658.01 

279.74 

M3-5 

F12 

863.01 

2225.94 

'  2672.40 

243.05 

♦ 

M4-1 

F13 

860.39 

2037.91 

'’^3101.05 

277.05 

* 

M4-2 

F14 

860.56 

3252.33 

* 

3040.17 

424.67 

M4-3 

F12 

847.90 

2919.14 

* 

3298.97 

454.17 

M4-4 

F13 

859.50 

2457.66 

* 

* 

F14 

859.94 

3065.20 

* 

♦ 

M5-1 

F14 

857.48 

1985.81 

2139.20 

237.86 

1748.98 

137.29 

M5-2 

F13 

853.51 

3095.07 

♦ 

3198.03 

357.07 

M5-3 

F13 

859.33 

2773.87 

3028.57 

314.05 

2768.52 

139.91 

M5-4 

F14 

857.49 

2219.99 

2266.20 

269.96 

2125.41 

127.56 

M6-1 

F14 

850.07 

2888.39 

♦ 

* 

M6-2 

F13 

854.04 

2551.80 

* 

* 

M6-3 

liBil 

846.97 

3158.72 

3131.14 

327.80 

3156.74 

M7-1 

13E1 

859.22 

2428.68 

2475.00 

308.03 

2364.49 

M7-2 

F12 

849.39 

2349.75 

2456.94 

300.83 

2357.10 

M7-3 

F13 

854.98 

3264.78 

* 

3269.50 

M8-1 

F14 

848.18 

2548.46 

'  3802.60 

431.60 

3380.70 

M8-2 

F13 

853.77 

2742.37 

* 

* 

M8-3 

F12 

845.71 

3014.86 

353.60 

3184.45 

M9-1 

F13 

852.12 

2766.80 

242.70 

2657.60 

MlO-1 

EEQ 

846.25 

3283.46 

3576.95 

373.36 

3291.50 

MlO-2 

1291 

841.57 

3446.84 

3637.09 

419.75 

3412.61 

189.27 

MlO-3 

F13 

857.97 

3010.32 

3209.01 

186.32 

3066.39 

97.81 

Ml  0-4 

F14 

853.77 

2305.00 

2376.51 

156.68 

2218.76 

70.71 

Ml  0-5 

F15 

842.59 

3105.89 

3166.29 

307.09 

3066.52 

138.30 

Ml  0-6 

F13 

857.48 

2801.96 

2868.91 

164.75 

2704.28 

89.00 

Ml  0-7 

F14 

853.74 

2325.14 

2437.31 

140.64 

2312.87 

74.31 

MlO-8 

F15 

846.36 

2161.87 

2170.73 

143.50 

2013.95 

66.07 

Cases 

35 

20 

26 

Avg.T,(K) 

2722.99 

2868.60 

±271.92 

2861.05 

±  217.46 

Relative  Error 

9.48% 

7.60% 

Avg.  ISR 
Uncertainty 

8.42% 

Notes:  *  N/A 

'  Based  upon  only  two  profiles 
^  SN  ratio  «  0.1 
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DMSP  vs.  ISR  Te  Comparison  over  Millstone  Hill 

Of  the  37  Millstone  conjunctions,  DMSP  Te  measurements  consistently  exceeded 
the  “filtered”  ISR  Te  estimates  for  all  three  pulse  lengths.  The  difference  between  DMSP 
Te  and  corresponding  ISR  Te  decreased  towards  solar  maximum.  Ten  values  deviated 
from  this  trend  (four  DMSP  and  six  ISR),  probably  due  to  SSIES  SENPOT  problems  and 
highly  uncertain  ISR  estimates  based  on  only  two  profiles  or  with  high  SN  ratios.  Figure 
25  below  shows  all  available  DMSP  and  ISR  Te  measurements  for  the  37  conjunctions. 
The  DMSP  spacecraft  involved  in  each  case  is  shown  as  the  last  term  of  the  case  number. 


DMSP  vs.  Millstone  Hill  ISR  Te  Comparison 
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Case  Number 


Figure  25.  DMSP  Te  vs.  All  Available  Millstone  Hill  ISR  Comparison  Values.  Note  that 
DMSP  Te  values  exceed  ISR  values  in  most  cases.  This  disparity  shrinks  as  POLITE 
cases  occur  closer  to  solar  maximum  (right  side  of  graph). 
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A  table  containing  the  DMSP  and  Millstone  Hill  ISR  Te  comparison  values  supporting 
Figure  25  can  be  found  in  Appendix  B. 

The  overall  average  percent  difference  between  the  DMSP  Te  data  and  all 
available  ISR  Te  estimates  for  33  of  the  37  cases  was  found  to  be  24.6  percent,  with  a 
standard  deviation  of  19.02  percent.  This  is  roughly  three  times  greater  than  the  average 
ISR  uncertainty  and  about  six  times  greater  than  the  DMSP  random  uncertainty  (3.71 
percent).  This  consistently  large  percent  difference  relative  to  ISR  and  DMSP 
uncertainties  could  suggest  a  bias  in  DMSP  or  ISR  Te  measurements.  DMSP  and  ISR  Te 
measurements  generally  become  closer  towards  solar  maximum,  with  percent  differences 
decreasing  from  about  46  percent  near  solar  minimum  to  15-20  percent  towards  solar 
maximum.  Table  13  below  shows  the  breakdown  of  the  average  DMSP  -  ISR  percent 
differences  and  standard  deviations  by  pulse  length.  Scatter  plots  and  tables  presenting 
each  case’s  percent  difference  by  pulse  length  are  found  in  Appendix  A. 


Table  13.  DMSP  -  Millstone  ISR  Average  Percent  Differences 


Pulse  Length 

#  of  Cases 

Avg  %  Diff. 

%  Diff.  a 

31 

18.95 

16 

14.12 

13.39 

2000  ps 

22 

19.39 

Total  Avg. 

24.60 

19.02 

The  1000  ps  cases  show  the  smallest  average  percent  difference  of  13.39  percent.  This 
result  can  be  attributed  to  the  fact  that  most  of  these  cases  occur  closer  to  solar  maximum, 
where  DMSP  Te  -  ISR  Te  agreement  improved. 


69 


Photoelectron  and  Solar  Cycle  Influences  on  DMSP  Te  Values 

Referring  back  to  Figure  25,  DMSP  Te  values  show  distinct  solar  zenith  angle 
(SZA)  and  solar  cycle  related  trends  that  are  independent  of  spacecraft.  Percent 
differences  are  smaller  for  conjunctions  with  a  larger  SZA,  including  those  cases 
occurring  diiring  darkness,  when  the  solar  zenith  angle  exceeded  118  degrees.  This  trend 
is  most  evident  during  POLITE  campaigns  1-5  in  particular,  as  DMSP  flyovers  with  a 
smaller  SZA  (corresponding  to  more  frequent  photoionization  and  more  photoelectrons 
produced)  resulted  in  much  higher  Te  values  than  both  corresponding  ISR  Te  values  and 
neighboring  DMSP  Te  values  within  the  same  POLITE  campaign.  Analysis  of  POLITE 
1-5  conjunctions  with  DMSP  Te  values  (from  Figure  25)  that  are  clearly  higher  than  the 
main  Te  grouping  revealed  an  average  SZA  of  66  degrees,  with  DMSP  Te  values 
averaging  46  percent  higher  than  average  ISR  Te  values  (See  Table  14).  On  the  other 
hand,  the  remaining  conjunctions  with  DMSP  Tg  values  closer  to  ISR  Te  values  during 
POLITEs  1-5  averaged  a  103-degree  SZA  and  19  percent  difference  (see  Table  15).  This 
trend  suggests  a  pronounced  photoelectron  impact  on  Te,  where  photoelectrons  are 
collected  by  the  electron  probe,  and  lead  to  elevated  Te  measurements.  This 
photoelectron  influence  on  DMSP  Te  measurements  decreased  towards  solar  maximum 
as  ambient  electron  densities  increasd  (see  discussion  on  pg.  33).  In  fact,  the  19  percent 
difference  from  those  cases  during  POLITEs  1-5  not  suspected  for  photoelectron 
contamination  was  quite  consistent  with  the  overall  percent  difference  results  (15-20 
percent)  for  all  cases  towards  solar  maximum. 
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Table  14.  Millstone  Hill  POLITE  1-5  Cases:  “Low”  Solar  Zenith  Angle 


Case  Number 

Satellite 

SZAH 

DMSP  Te 

ISR  Avg.  Te 

DMSP  -  ISR 

M2-1 

66 

4988.4 

2958.0 

2030.4 

M2-2 

F12 

65 

4522.9 

2600.3 

1922.6 

M2-3 

F12 

65 

4716.7 

3181.5 

1535.2 

M3-1 

F12 

32 

2997.1 

1359.6 

M3-3 

71 

2530.6 

F13 

67 

3856.4 

2658.0 

1198.4 

M4-2 

F12 

75 

4223.8 

3146.3 

M4-3 

F14 

70 

4146.2 

M4-5 

F14 

75 

M5-2 

F13 

71 

4084.4 

3146.6 

937.9 

M5-3 

F13 

69 

2857.0 

964.7 

Mean 

66 

2931.8 

1358.5 

Std.  Deviation 

11.8 

355.5 

236.7 

370.7 

%  Difference 

46.34% 

Table  15.  Millstone  Hill  POLITE  1-5  Cases:  “High”  Solar  Zenith  Angle 


Case  Number 

Satellite 

szA  n 

DMSP  Te 

ISR  Avg.  Te 

DMSP  -  ISR 

*M1-1 

F12 

2921.0 

713.1 

Ml-2 

F12 

98 

2809.7 

mEMM 

-301.2 

M3-2 

wsm 

102 

2989.9 

221.6 

M3-5 

msm 

102 

mmm 

2449.2 

963.6 

M4-1 

F13 

99 

2960.5 

2569.5 

391.0 

M4-4 

F13 

99 

2457.7 

445.4 

M5-1 

F14 

97 

1956.0 

746.6 

M5-4 

F14 

97 

2822.4 

2203.9 

618.5 

Mean 

102.8 

2968.0 

2493.1 

474.8 

Std.  Deviation 

10.4 

233.1 

394.8 

390.0 

%  Difference 

19.05% 

Note:  '  Spacecraft  in  Darkness 


DMSP  vs.  ISR  Te  Comparison  over  Sondrestrom 

Comparison  of  the  six  Sondrestrom  conjunctions  shows  a  similar  trend  of  DMSP 
Te  values  exceeding  its  corresponding  extrapolated  ISR  450  ps  Te  value  in  all  but  one 
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instance.  The  overall  DMSP  Te  percent  difference  from  ISR  values  was  19.15  percent 
with  a  standard  deviation  of  25.23  percent  (see  Table  16).  This  smaller  percent 


Table  16.  DMSP  -  ISR  Te  Comparison  over  Sondrestrom 


Case# 

Satellite 

Solar  Zenith 
Angle 

DMSP  Te  (K) 

a 

450ps  Ext.  Te  (K) 

%  Diff. 

‘S2-1 

F12 

2160.0 

361.77 

N/A 

N/A 

S3-1 

F14 

76° 

237.12 

3134.94 

28.76 

S3-2 

63° 

332.57 

3527.03 

22.46 

S5-1 

F14 

490 

H||||||||||EQQ|^| 

122.80 

3359.83 

13.82 

^S6-l 

mm 

118° 

2254.97 

-21.06 

S8-1 

mm 

103° 

2726,7 

1740.20 

56.69 

S9-1 

F13 

103^ 

504.65 

1838.87 

14.20 

Avg.  %  Diff. 

19.15 

a  of  Avg.  %  Diff. 

25.23 

Notes:  '  jH  =  Spacecraft  in  darkness 

^  Terminator  Zenith  Angle  =118° 


difference  compared  to  the  average  Millstone  Hill  percent  difference  of  24.6  percent 
could  be  misleading  due  to  the  small  number  of  Sondrestrom  cases  involved  and  the  large 
uncertainty  involved  with  extrapolating  Sondrestrom  ISR  Te  data.  Both  DMSP  and 
corresponding  ISR  Te  measurements  are  higher  in  cases  during  summer  in  daylight  with 
smaller  solar  zenith  angles,  than  in  the  cases  during  winter  with  larger  solar  zenith  angles. 
Figure  26  shows  a  graphical  representation  of  these  conjunctions.  Note  that  the  DMSP 
spacecraft  involved  in  each  case  is  shown  as  the  last  term  of  the  case  number.  In  Case 
6-1,  the  one  instance  where  the  DMSP  Te  is  less  than  the  extrapolated  ISR  Te  value,  the 
altered  mode  of  SSDES  measurement  collection  at  the  time  (POLITE  6  -  see  pg.  38,  Item 
2)  could  have  contributed  to  the  overly  smoothed  SSBES  Tg  values.  The  fact  that  the 
SSIES  and  ISR  Tg  values  trace  similar  Tg  patterns  for  five  of  the  six  conjunctions 
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suggests  that  both  instruments  are  equally  responsive  to  the  surrounding  ionospheric 
phenomena,  and  perhaps  a  similar  DMSP  Te  measurement  bias  exists  at  high  latitudes  as 
well.  Due  to  only  six  bona-fide  Te  comparisons  examined  over  Sondrestrom,  and  the 
expected  large  variability  of  high  latitude  Te  over  short  distances,  it  is  difficult  to  make 
any  further  assessments  without  pursuing  a  more  comprehensive  study. 


DMSP  vs.  Sondrestrom  ISR  Te  Comparison 


Figure  26.  Conjunctions  of  DSMP  -  450  ps  ISR  Extrapolated  Te  Values  Above 
Sondrestrom.  DMSP  values  are  averaged  within  a  one-degree  circle  “up  B”  centered  at 
65.9  N,  31 1.2  E  (Geographic).  Note  that  ISR  extrapolation  was  not  possible  for  Case  2-1. 


Unusual  DMSP  Te  Profiles 

Closer  inspection  of  DMSP  Te  behavior  along  the  44  orbits  containing  the 
Millstone  Hill  and  Sondrestrom  POLITE  conjunctions  reveal  possible  instrument  and 
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data  processing  related  anomalies  that  could  have  distorted  averaged  Te  values  from  14  of 
the  conjunctions.  First,  Dr.  Hairston  cautioned  that  Te  values  become  unreliable  when 
the  SSIES  SENPOT  cannot  remain  at  the  plasma  potential,  and  consequently  drops  to  the 
spacecraft  potential  [Hairston,  private  correspondence,  21  Dec  00] .  This  seemed  to  be 
the  case  most  noticeably  during  a  considerable  sector  of  the  daylight  portion  of  the  FI  4 
orbits  during  POLITE  10,  and  more  in  general  from  where  the  solar  zenith  angle 
decreases  from  its  maximum  to  minimum  value,  from  about  140  to  40  degrees.  As  soon 
as  the  spacecraft  crossed  the  terminator  into  daylight,  Te  values  dropped  to  between  700  - 
1000  K.  Electron  temperatures  then  abruptly  jumped  to  reasonable  once  again  when  the 
solar  zenith  angle  began  climbing  from  its  minimum  value  (~  40  degrees).  Once  the 
zenith  angle  reached  its  maximum  value  (~  140  degrees)  and  began  decreasing,  Te  values 
gradually  decreased  seemingly  more  in  response  to  an  impredicted  SENPOT  increase 
than  to  true  ionospheric  variations.  Figure  27  shows  this  pattern  of  Te  behavior  for  a  F14 
orbit  during  POLITE  10  that  contains  Millstone  Hill  case  Ml  0-1.  F12  SENPOT 
problems  seem  to  impact  two  Millstone  conjimctions  during  earlier  POLITE  campaigns 
as  well,  but  not  as  dramatically  as  with  the  F14  POLITE  10  cases  (see  Chapter  m. 

Table  8). 
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DMSP  Orbit  -  POLITE  10  (F14) 
1  JulOO  1310-  1452Z 
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Figure  27.  F14  Orbit  Containing  Millstone  Hill  Conjunction  MlO-1.  As  the  spacecraft 
traverses  the  sunlit  portion  of  its  orbit  (much  of  the  Northern  Hemisphere),  SENPOT 
goes  missing  and  Te  values  drop  to  nominal  values  of  700  -  lOOOK.  Average  DMSP  Te 
forMlO-1  is  81 1.0  K. 

Playing  more  of  a  minor  impact  on  results,  as  mentioned  in  Chapter  HI,  the  SSIES 
electron  probe  was  in  a  retarded  Te  measurement  collection  mode  during  the  entire 
POLITE  6  time  frame  (22  -  25  November,  1998)  [Hairston,  private  correspondence,  21 
Dec  00].  This  reduced  the  resolution  of  the  Te  structure  and  resulted  in  a  smoother,  step¬ 
like  Te  structure  through  the  entire  spacecraft  orbit  (see  Figure  28).  As  another  possible 
cause  for  skewed  results.  Dr.  Sultan  of  AFRL  advised  us  to  be  wary  of  F13  Te 
measurements  being  too  high  during  the  evening  pass  over  Millstone  Hill  [Sultan,  private 
correspondence,  26  Sep  00].  Table  17  shows  a  breakdown  of  conjunctions  possibly 
affected  by  instrument  or  data  processing  related  anomalies.  Note  that  the  four  cases 
(M8-3,  MlO-1,  MlO-4,  MlO-7)  that  show  a  SENPOT  anomaly  as  a  significant  impact  on 
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the  resulting  DMSP  mean  Te  were  excluded  from  our  results  statistics.  A  significant 
impact  is  defined  as  an  anomaly  producing  a  DMSP  Tg  value  more  than  1500  K  away 
from  corresponding  ISR  Tg  estimates.  Evidence  of  anomaly  impact  on  a  case’s  DMSP 
electron  temperature  profile  is  best  viewed  in  the  DMSP  Tg  plots  for  all  37  Millstone 
conjunction  flyovers,  which  are  available  as  a  separate  document  upon  request. 


DMSP  Orbit  -  POLITE  6  (F12) 
23  Nov  98  1335-  1517Z 


6000 
5000 
4000 

S' 

^  3000 

o 

*”  2000 
1000 
0 

000)CMT-c000l^'t00C0T-0)lO0)'^C000O 
CO  lO  CD  CO  T-  ' 

Geographic  Latitude 


Figure  28.  DMSP  Tg  Measurements  from  an  Orbit  During  POLITE  6.  A  smoother, 
step-like  profile  results  from  a  slower  measurement  collection  mode. 
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Table  17.  Cases  with  Possible  Distorted  Te  Results 


Case 

Satellite 

Cause  of  Anomaly 

Impact  on  Te 

Impact  on  ISR 
Comparison 

M3-2 

F12 

SENPOT 

Sporadic  Fluctuations 

Moderate 

^^9 

F12 

SENPOT 

Sporadic  Fluctuations 

Significant 

DBffH 

F14 

SENPOT 

Te  700- 1000  K 

F14 

SENPOT 

Te  700- 1000  K 

F14 

SENPOT 

Te  700- 1000  K 

SSIES  settings 

Te  too  high 

Moderate 

SSIES  settings 

Te  too  high 

Moderate 

WHHI 

F13 

SSIES  settings 

Te  too  high 

Moderate 

F13 

SSIES  settings 

Te  too  high 

Moderate 

F13 

SSIES  settings 

Te  too  high 

Moderate 

F14 

Slower  Te  Collection  Mode 

Smoother  Profile 

Minor 

F13 

Slower  Te  Collection  Mode 

Smoother  Profile 

Minor 

F12 

Slower  Te  Collection  Mode 

Smoother  Profile 

Minor 

S6-1 

F14 

Slower  Te  Collection  Mode 

Smoother  Profile 

Minor 

Note:  ’  Cases  Discarded  for  Unreliable  Te  Values. 


DMSP  Te  Random  Noise  Level  Results 

After  four  of  the  DMSP  Te  profiles  were  removed  for  probable  SENPOT  error, 
the  average  noise  level,  calculated  as  the  standard  deviation  as  a  fraction  of  the  mean 
(SDFM),  of  the  remaining  33  data  sets,  turned  out  to  be  3.71  percent.  This  result  falls 
well  within  the  quoted  SSIES  Te  measurement  accuracy  of  10  percent  in  Rich’s  SSIES 
User’s  Guide  [1994]  despite  the  rudimentary  assumption  of  a  linear  behavior  of  mid 
latitude  electron  temperatures  to  de-trend  the  data.  Table  18  contains  a  summary  of 

SDFM  results  by  case,  where  Te,  =  Filtered  Mean  Te  and  CTr\  =  Standard  Deviation  of 
Filtered  Residuals. 
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Table  18.  DMSP  Random  Noise  Level  (SDFM)  Calculations 


Case  # 

Satellite 

Solar  Zen. 
Angle 

Tot. 

Pts. 

Pts. 

Removed 

Te,  (K) 

crr^ 

SDFM 

liiiSM 

F12 

1 

2895.0 

113.233 

3.91% 

F13 

98 

29 

2 

145.950 

5.21% 

F12 

66 

31 

1 

5076.0 

F12 

65 

41 

1 

4535.5 

65 

36 

0 

4716.7 

32 

18 

2 

4367.5 

158.800 

3.64% 

F12 

13 

1 

3345.8 

286.755 

M3-3 

F13 

71 

29 

2 

4097.0 

122.360 

2.99% 

msM 

F13 

67 

22 

1 

3833.3 

121.622 

3.17% 

ms^ 

102 

39 

1 

3425.8 

189.208 

5.52% 

M4-1 

F13 

99 

40 

1 

2952.3 

124.063 

4.20% 

F12 

75 

26 

0 

4223.8 

172.367 

F14 

39 

2 

4175.1 

150.046 

■ee 

M4-4 

F13 

99 

26 

1 

2912.0 

87.275 

3.00% 

■EH 

75 

38 

1 

4421.1 

153.842 

■BE 

BHI 

97 

39 

1 

2709.5 

107.161 

■ee 

F13 

71 

27 

2 

4092.0 

57.119 

1.40% 

F13 

69 

23 

ESI 

F14 

97 

34 

1 

F14 

11 

1 

2976.0 

B3EB 

97 

40 

3 

3503.3 

52.422 

■BEI 

msm 

70 

37 

0 

3693.5 

99.322 

F14 

41 

3 

2570.5 

232.428 

9.04% 

msm 

27 

IIESEI 

30.061 

1.16% 

msam 

86 

24 

1 

357.924 

10.09% 

■IgMB 

F14 

15 

0 

3129.3 

66.524 

2.13% 

M8-2 

F13 

99 

24 

1 

97.064 

3.25% 

'M8-3 

F12 

75 

22 

2045.891 

36.58% 

M9-1 

F13 

98 

28 

2 

77.985 

2.34% 

UBiMl 

F14 

40 

29 

1 

805.7 

66.780 

8.29% 

F15 

33 

39 

1 

4147.9 

91.625 

2.21% 

MlO-3 

F13 

68 

1 

102.686 

2.96% 

F14 

95 

^El 

1 

771.4 

33 

12 

0 

M10--6 

F13 

68 

mm 

1 

90.960 

2.74% 

'MlO-7 

F14 

95 

E^l 

1 

768.7 

59.210 

7.70% 

MlO-8 

F15 

101 

IEH 

1 

63.944 

2.25% 

Avg. 

SDFM 

3.71% 

Notes:  *  Cases  not  used  in  SDFM  average  due  to  SENPOT  error 


Spacecraft  crosses  terminator  solar  zenith  angle  (118°) 
^  Probable  sporadic  SENPOT  error 
^  =  Spacecraft  in  dark 
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As  the  table  shows,  most  individual  noise  calculations  are  under  five  percent.  Six 
of  the  nine  SDFM  values  over  five  percent  can  be  at  least  partially  attributed  to  SENPOT 
problems.  No  distinct  trend  could  be  detected  with  respect  to  season,  however  noise 
levels  are  generally  about  1-2  percent  smaller  closer  to  solar  max  and  for  those  cases 
where  the  spacecraft  is  in  the  dark. 

Although  a  linear  fit  seemed  reasonable  to  de-trend  most  data  sets,  as  in  the 
example  in  Chapter  HI,  a  few  cases  showed  a  variable  Te  structure  where  a  straight  line 
fit  was  not  appropriate.  Also  the  2ar  threshold  established  to  remove  outlying  points  did 
not  always  do  so  after  just  one  de-trending  and  filtering  iteration.  In  fact,  a  poor  linear  fit 
after  one  iteration  of  filtering  highly  variable  data  can  still  contain  points  which  should  be 
outliers  and  can  also  exclude  points  that  are  part  of  the  physical  trend.  This  can  result  in 
noise  estimates  that  are  not  representative  of  the  true  data  variation.  Figures  29  and  30 
shows  the  linear  fits  of  raw  and  filtered  Te  data  for  Case  M7-1  that  demonstrate  these 
shortcomings.  Further  inspection  of  this  case  suggests  that  the  SSDES  SENPOT  abruptly 
caused  unreliable  DMSPT  Te  values  when  the  spacecraft  crossed  into  the  sunlit  sector  of 
its  orbit  near  46  degrees  N  (geographic).  Thus  a  subjective  “human”  analysis  for  such 
cases  would  help  complement  the  aforementioned  noise  level  determination 
methodology. 
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Figure  29.  Initial  Fit  of  Case  M7-1  DMSP  Te  Data  Showing  Three  Outliers. 


Figure  30.  Poor  Linear  Fit  of  “Filtered”  Data.  Three  of  the  nine  points  outside  of  the 
recalculated  2<jr  boundary  seem  to  be  “bad”  data,  whereas  the  other  six  seem  to  follow 
the  physical  trend  of  the  data. 
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In  instances  such  as  these,  a  simple  curve  fit  and  a  second  filtering  iteration  could  have 
kept  and  removed  the  appropriate  points.  Furthermore,  applying  Chauvenet’s  Criterion 
to  reject  unreasonable  measurements  might  have  more  accurately  removed  the 
appropriate  outliers  than  our  established  lor  boundary. 

Case  Comparisons 

Two  case  comparisons  were  performed  to  evaluate  the  consistency  of  DMSP  Te 
measurements  collected  by  different  spacecraft  passing  over  Millstone  Hill  within 
minutes  of  each  other.  Both  comparisons  involve  F14  and  F12  spacecraft  and  are  listed 
in  Table  19. 


Table  19.  Case  Comparison  Details 


Sat. 

Pts.  w/in  5° 

F14 

860.56 

75° 

39 

4146.2 

Compared  to: 

B^^B 

F12 

847.90 

o 

O 

26 

4223.8 

99282 

0059 

F14 

859.22 

' 

41 

2486.8 

Compared  to: 

M7-2' 

99282 

0114 

F12 

849.39 

27 

2580.7 

Notes:  *  Bf  =  Spacecraft  in  dark 

Crosses  terminator  zenith  angle  (118°)  within  five-degree  circle 


In  both  instances,  FI 2  flies  over  Millstone  Hill  about  20  minutes  later  at  a  slightly  lower 
altitude  than  F14  and  yields  Te  values  about  80  -  90  K  higher  than  those  from  F14.  This 
Te  difference  is  within  the  random  noise  level  (3.71  percent)  calculated  during  this 
research.  In  addition,  we  can  probably  assume  steady-state  ionospheric  conditions  during 


81 


these  two  periods.  Therefore,  the  SSIES  electron  probes  on  both  spacecraft  seem  to  be 
taking  consistent  Te  measurements  during  these  periods.  Figure  31  shows  how  cases 
M7-1  and  M7-2  trace  similar  Te  profiles  approaching  the  terminator  firom  darkness  into 
sunlight.  Note  that  Te  data  was  missing  past  47°  N  (Geographic)  for  Case  M7-2. 


F14  DMSP  Te  ±3nnin  of  Millstone  Hill 


«  ......  Avg.  Te  within  5°=  2486.8K 

Geog.  Latitude  388,73k 


F12  DMSP  Te  ±  3min  of  Millstone  Hill 


Figure  31.  DMSP  Te  Profile  Comparison  of  (a)  Case  M7-1  against  (b)  Case  M7-2.  The 
mean  DMSP  Te  for  Case  M7-2  within  the  five-degree  circle  exceeds  that  fi'om  Case  M7-1 
by  only  93.9  K. 
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Case  Studies  -  Mid-Latitude  Ionospheric  Phenomena 

Evidence  of  characteristic  mid-latitude  ionospheric  features  such  as  subauroral 
ion  drifts  (SAID),  electron  density  (ne)  troughs  (MT),  and  light  ion  troughs  (LIT)  is 
present  in  several  Millstone  conjunction  profiles.  After  inspecting  a  20-degree  latitude 
range  of  SSIES  data  surrounding  Millstone  Hill  during  the  37  conjunctions.  Table  20 
below  gives  a  summary  of  the  features  that  occurred.  Note  that  bad  or  missing  RPA  and 
SENPOT  data  limited  the  completeness  of  this  effort. 


Table  20.  Case  Summary  of  Mid  Latitude  Features 


Case 

SAID 

MT 

LIT 

Case 

SAID 

MT 

LIT 

Ml-1 

X 

X 

X 

M6-1 

X 

HiB 

X 

X 

X 

X 

msm 

M2-2 

X 

M7-1 

M2-3 

M7-2 

M3-1 

M7-3 

X 

EB 

mmm 

eb 

X 

1  M3-4 

eb 

M9-1 

HEB 

X 

X 

Ml  0-1 

X 

X 

laiiiii 

X 

■laiiiHM 

X 

MW 

X 

X 

Ml  0-4 

X 

BEB 

X 

laifaai 

X 

msm 

X 

MlO-7 

MlO-8 

B£B 

Being  located  equatorward  of  the  auroral  oval,  Millstone  Hill’s  location  was  ideal 
to  detection  of  at  least  part  of  a  mid-latitude  Ue  trough.  Mid  latitude  troughs  were  mainly 
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identified  by  correlating  a  sloped  ne  profile  with  a  corresponding  inversely  related  Te 
profile.  Figure  32  shows  an  example  of  a  MX  identified  during  a  DMSP  flyover  during 
POLITE  10  (Case  MlO-3).  Electron  temperature  seems  to  reach  a  general  maximum 
where  the  Ue  trough  bottoms  out. 


Figure  32.  Case  Ml  0-3,  DMSP  Plots  of  (a)  Te  and  (b)  ne  Showing  Support  for  a  Mid- 
Latitude  Trough  around  46°N  Geographic  Latitude. 


The  best  indication  of  a  SAID  and  LIT  were  shown  in  Case  Ml-1 .  Recall  that  a 


SAID  is  typically  a  region  of  a  rapid  westward  ion  drift  which,  via  upward  thermal 
conduction,  can  cause  elevated  Te  values  at  DMSP  altitude.  Meanwhile,  recollect  that  the 
LIT  is  rapid  density  depletion  of  kT  or  He^  that  is  most  pronounced  in  the  nighttime 
winter  evening  sector.  All  of  these  ingredients  exist  to  support  a  SAID,  LIT,  and  MT 
occurrence  during  Case  Ml-1  (See  Figure  33).  Note  that  in  the  ion  velocity  plot  (b),  +Vy 
is  oriented  to  the  left  of  the  spacecraft’s  forward  direction  of  motion,  which  for  this  case 
corresponds  to  a  westward  drift  [Hairston,  private  correspondence,  23  Aug  00] . 
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F12  DMSP  Ion  Density  +  3min  of  Millstone  Hll 


Geog.  Latitude 


F12  DMSP  Ion  Velocity  Components  ±  3min  of  Millstone  Hill 


a. 


F1 2  DMSP  Te  ±  3min  of  Millstone  Hill 


F12  DMSP  Electron  Density  ±  3min  of  Millstone  Hill 


C. 


b. 


d. 


Figiire  33.  Simultaneous  SSEES  plots  of  (a)  Ion  Density,  (b)  Ion  Velocity,  (c)  Te,  and  (d) 
ne  profiles  during  Case  Ml -1  on  14  Feb  96  Showing  Evidence  of  a  LIT,  SAID,  and  MT 
towards  the  Right  Side  of  the  Plots. 
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V.  Summary.  Conclusions  and  Recommendations 


Summary 

Geomagnetic  disturbances  in  the  near  earth  space  environment  can  adversely 
affect  numerous  military  and  DOD  systems  and  operations.  To  improve  the  prediction 
accuracy  of  such  disturbances,  the  next  generation  of  space  environment  forecast  models 
(i.e.  GAIM)  aims  to  incorporate  near  real-time  space-based  and  ground-based 
measurements  into  a  physically  realistic  model  field.  DMSP  SSIES  electron 
temperatures  represent  one  such  space-based  parameter  to  be  ingested  into  GAIM.  Thus 
the  need  arose  to  validate  and  assess  the  reliability  and  noise  level  of  DMSP  SSIES 
electron  temperatures. 

SSIES  Te  measurements  were  compared  against  near  simultaneous  Millstone  Hill 
and  Sondrestrom  ISR  Te  measurements  taken  during  POLITE  Campaigns  1-10  spanning 
from  Winter  1996  to  Summer  2000.  Incoherent  Scatter  Radar  Te  measurements  were 
filtered  to  remove  highly  erroneous  or  noisy  data  and  then  either  averaged  or  extrapolated 
to  DMSP  altitude  to  determine  a  comparison  value.  Conjunction  criteria  similar  to  that 
used  in  the  Sultan  and  Rich  [2000]  recent  ion  density  study  led  to  37  Millstone  Hill 
comparison  cases,  while  a  more  restrictive  criteria  resulted  in  seven  possible  Sondrestrom 
cases.  Percent  differences  were  then  calculated  between  each  case’s  Te  comparison 
values. 

A  first  guess  linear  regression  was  then  performed  on  the  37  DMSP  Millstone  Hill 
flyovers  to  estimate  random  noise  level  as  the  standard  deviation  as  a  fraction  of  the 
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mean.  This  approach  was  not  appropriate  to  approximate  the  highly  variable  Tg  structure 
at  auroral  latitudes,  thus  a  noise  level  was  not  computed  for  the  Sondrestrom  flyovers. 

Conclusions 

For  both  locations,  and  at  all  available  ISR  pulse  lengths,  DMSP  Tg  values 
consistently  exceeded  comparable  averaged  ISR  values.  For  the  37  Millstone  Hill 
conjunctions,  DMSP  Tg  values  averaged  24.6  percent  greater  than  Millstone  ISR  values. 
This  result  is  nearly  three  times  greater  than  the  average  computed  Millstone  ISR 
uncertainty.  Percent  differences  generally  were  largest  when  the  solar  zenith  angle  was 
smallest,  especially  during  solar  minimum  when  ambient  electron  densities  were  low, 
suggesting  a  pronounced  photoelectron  impact  on  DMSP  Tg  values.  Satellite  versus  radar 
Tg  percent  differences  for  the  Millstone  Hill  conjunctions  decreased  as  the  solar  cycle 
approached  solar  maximum.  Also,  percent  differences  were  generally  smaller  when  the 
conjunction  region  was  in  darkness;  perhaps  due  to  photoelectrons  collected  by  the 
SSIES  electron  probe  (when  svmlit)  interfering  with  the  Tg  measurement  procedure  and 
skewing  resulting  Tg  values. 

For  the  six  Sondrestrom  conjunctions,  DMSP  Tg  values  averaged  19  percent 
greater  than  ISR  measurements.  Although  this  average  percent  difference  is  smaller  than 
that  calculated  for  the  Millstone  Hill  eonjunctions,  the  Sondrestrom  percent  difference 
tally  showed  a  four  percent  larger  standard  deviation.  The  smaller  average  percent 
difference  result  over  Sondrestrom  must  be  viewed  in  light  of  several  faetors.  These 
include  the  highly  dynamic  nature  of  the  high  latitude  ionosphere,  the  ISR  Tg  filtering  and 
extrapolation  routine  involved  working  with  large  errors  and  uncertainties,  and  the  fact 
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that  only  seven  conjunctions  were  considered.  Consequently,  a  more  thorough  Te 
comparison  for  Sondrestrom  would  likely  generate  more  trustworthy  statistics. 
Meanwhile,  the  concurrent  electron  density  study  over  the  two  ISR  sites  revealed  average 
DMSP  lie  measurements  that  were  lower  than  average  ISR  ne  values,  thus  verifying  the 
expected  inverse  Te  vs.  ne  relationship  for  a  given  amount  of  energy  [Keyser,  Thesis 
Draft.  2001]. 

Average  DMSP  random  noise  over  Millstone  Hill  for  the  37  cases  was  computed 
at  3.71  percent,  which  is  well  within  the  published  SSEES  Te  measurement  uncertainty  of 
ten  percent.  [Rich,  1994],  Individual  conjunction  noise  estimates  that  were  larger  than 
six  percent  were  most  likely  due  to  a  SENPOT  problem  disrupting  the  Te  profile,  or  when 
our  assumption  of  a  smooth,  linear  Te  variation  over  Millstone  Hill  was  invalid. 

Major  sources  of  error  for  this  study  included  instrument  anomalies  existing 
during  parts  of  DMSP  orbits,  the  random  DMSP  noise  level,  the  uncertainty  involved 
with  averaging  or  extrapolating  ISR  values  to  DMSP  altitudes,  and  the  inherent  Te 
uncertainty  already  accompanying  the  ISR  Te  values,  especially  at  higher  altitudes.  The 
fact  that  the  ISR  Te  values  relied  upon  a  number  of  key  assumptions,  together  with  the 
frequent  highly  uncertain  ISR  Te  behavior  at  DMSP  altitude,  showed  that  such 
measurements  were  not  quite  as  reliable  a  standard  as  anticipated.  We  accounted  for  the 
aforementioned  sources  of  error  as  best  we  could  by  filtering  out  “bad”  ISR  Te  values  and 
removing  four  Millstone  conjunctions  from  our  calculations  due  to  clearly  unreliable 
DMSP  Te  values.  Thus,  we  remain  confident  in  our  findings.  Given  that  SSIES 
measured  Te  consistently  exceed  near-simultaneous  ISR  measured  Te,  and  that  different 
spacecraft  with  similar  flyover  times  and  solar  zenith  angles  over  Millstone  Hill  show 
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comparable  Te  values  (within  the  computed  random  DMSP  noise  level),  there  is  a 
possibility  that  some  sort  of  DMSP  Te  or  ISR  Te  bias  exists. 

Recommendations 

This  research  has  only  scratched  the  surface  of  accurately  validating  SSEES 
measured  parameters  to  be  ingested  into  GAIM.  The  first  step  fi’om  here  is  to  coordinate 
our  existing  results  with  experts  at  UTD,  AFRL,  and  Millstone  Hill  to  discover  the  source 
of  the  apparent  Te  bias  between  ISR  and  the  DMSP  SSIES.  This  will  surely  involve 
further  ISR  vs.  SSIES  data  comparisons  encompassing  other  latitude  ranges  and  a  wider 
range  of  geomagnetic  activity.  Incoherent  Scatter  Radar  data  from  Arecibo,  Puerto  Rico 
and  Jicamarca,  Peru  would  make  low  latitude  comparisons  possible,  while  Eiuropean 
Incoherent  Scatter  (EISCAT)  data  from  northern  Scandinavian  stations,  and  a  more 
thorough  Sondrestrom  SSIES  -  ISR  comparison,  would  expand  the  comparison  to  high 
latitudes.  Since  most  cases  in  this  study  occurred  during  geomagnetic  quiet  conditions 
(Kp  <  3),  examining  cases  during  periods  of  more  active  geomagnetic  storm  levels  will 
broaden  the  comparison  database  as  well.  With  a  more  comprehensive  database  of 
DMSP  vs.  ISR  measurement  comparison  results  with  respect  to  latitude,  solar  cycle,  solar 
illumination,  and  storm  level,  any  bias  that  might  exist  should  emerge. 

Since  ISR  returns  become  quite  weak  at  higher  altitudes,  especially  those  returns 
above  700  km  and  with  SN  ratios  less  than  0.1,  estimating  an  ISR  Te  value  at  DMSP 
altitude  involved  a  large  uncertainty.  Consequently,  more  in-situ  measurements  by 
rockets  or  other  spacecraft  need  to  be  performed  to  reduce  measurement  uncertainty  at 
DMSP  altitude. 
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Our  DMSP  Te  noise  level  study  was  limited  to  a  linear  fit,  mid-latitude  estimate. 
Ideally,  random  noise  error  needs  to  be  approximated  with  respect  to  latitude,  solar 
illumination,  solar  cycle  and  geomagnetic  activity.  A  start  toward  this  goal  would  be  to 
de-trend  and  compute  noise  levels  from  Te  data  from  other  latitude  sectors  of  the  1 147 
DMSP  orbits  available  for  this  study.  As  a  start  in  that  direction,  upon  sampling  Te  data 
fi'om  various  sections  of  DMSP  orbits,  the  equation:  y  =  a  +  bx^  best  fit  the  data  in 
nearly  one  third  of  32  cases.  GAIM  best  models  physical  parameters  exhibiting  slowly 
varying  ionospheric  behavior  as  at  low  and  mid-latitudes.  At  high  latitudes,  the  Kalman 
filter  would  treat  highly  varying  ionospheric  parameters  over  short  distances  as  noise 
[Sojka,  private  correspondence,  15  Nov  00] .  Thus  a  more  spatially  sensitive  routine  must 
be  considered  to  accoimt  for  the  short  scale  lengths  of  high  latitude  ionospheric 
phenomena. 

With  the  trend  of  next  generation  space  environment  models  relying  upon 
assimilating  spaced  based  measurements,  it  is  imperative  to  encourage  future  SSIES 
parameter  validations  and  to  continue  to  promote  near-real  time  telemetry  of  such 
parameters.  A  better  understanding  of  SSIES  characteristics  will  greatly  benefit  AFRL 
modeling  efforts  as  well  as  future  research  involving  DMSP  space  environment  data. 
Since  GAIM  has  the  capacity  to  ingest  new  data  types  as  they  become  available, 
expanding  the  stock  of  ionospheric  measurements  can  only  benefit  the  GAIM  initiative. 
As  a  result,  future  space  environment  forecasts  models  such  as  GAIM  will  become  more 
accurate  at  predicting  impending  geomagnetic  disturbances  that  could  impact  military 
and  DoD  systems  and  operations. 
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Case# 

Sat. 

DMSP  Mean 
Te(K) 

410  |xs 

Te(K) 

mmSt 

mmSt 

Ml-1 

F12 

2921.0 

2207.91 

32.30 

No  Data 

3602.15 

Ml-2 

F13 

2809.7 

* 

No  Data 

3110.90 

■sai 

F12 

4988.4 

2957.98 

* 

* 

F12 

4522.9 

2600.31 

♦ 

* 

M2-3 

F12 

4716.7 

2890.01 

63.21 

* 

3473.57 

M3-1 

4356.7 

3088.92 

* 

M3-2 

W3IMi 

3345.8 

2905.32 

* 

M3-3 

F13 

4071.7 

2530.55 

* 

* 

F13 

3856.4 

* 

* 

F12 

3412.8 

2225.94 

53.32 

2672.40 

27.71 

* 

M4-1 

2037.91 

45.27 

3101.05 

-4.53 

BBIKi 

4146.2 

^Em 

^  * 

F14 

4223.8 

44.69 

M4-4 

F13 

2903.1 

2457.66 

18.12 

* 

♦ 

F14 

4404.2 

3065.20 

43.68 

* 

* 

F14 

2702.6 

1985.81 

36.10 

2139.20 

26.34 

1748.98 

F13 

4084,4 

3095.07 

31.96 

* 

3198.03 

F13 

3821.7 

2773.87 

37.78 

3028.57 

26.19 

2768.52 

M5-4 

F14 

2822.4 

2219.99 

27.14 

2266.20 

24.54 

2125.41 

M6-1 

F14 

2974.5 

2888.39 

2.98 

* 

* 

F13 

3455.0 

2551.80 

35.39 

* 

* 

mm 

F12 

3693.5 

3158.72 

16.93 

3131.14 

msE 

MIA 

2486.8 

2428.68 

2.39 

2475,00 

2364.49 

ms^ 

F12 

2580.7 

2349.75 

9.83 

2456.94 

5.04 

2357.10 

■39 

3505.8 

3264.78 

7.38 

* 

3269.50 

M84 

F14 

3129.3 

2548.46 

22.79 

3802.60 

-17.71 

3380.70 

M8^2 

F13 

3005.8 

2742.37 

9.61 

* 

* 

‘M8-3 

F12 

5592.7 

3205.14 

3014.86 

3184.45 

M9-1 

F13 

3318.6 

2511.42 

32.14 

2766.80 

19.94 

2657.60 

‘MlO-l 

F14 

811.0 

3283.46 

3576.95 

3291.50 

Ml  0-2 

F15 

4144.6 

3446.84 

20.24 

3637.09 

13.95 

3412.61 

MlO-3 

F13 

3469.0 

3010.32 

15.24 

3209.01 

8.10 

3066.39 

'MlO-4 

F14 

779.1 

2305.00 

2376.51 

2218.76 

MlO-5 

F15 

3895.0 

3105.89 

25.41 

3166.29 

23.01 

3066.52 

MlO-6 

F13 

2801.96 

2868.91 

15.22 

2704.28 

*M10-7 

F14 

772.0 

2325.14 

2437.31 

2312.87 

MlO-8 

F15 

2839.4 

2161.87 

31.34 

2170.73 

30.80 

2013.95 

Cases 

32 

17 

Avg.  % 
Diflf. 

32.65 

18.32 

Avg.  % 
Difif.a 

20.15 

21.63 

Notes;  *  Te  Comparison  Value  Not  Available 

'  Percent  Difference  Not  Calculated  Due  to  Suspected  SENPOT  Error 
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Percent  Difference:  DMSP  Te  vs.  Millstone  Hill  410  |iS  ISR  Te 
(31  of  37  POLITE  Conjunctions  Represented) 


Case  Number 


Figure  34.  DMSP  vs.  Millstone  Hill  410  |j,s  ISR  Te  Percent  Difference  Case  Summary 


Figure  35.  DMSP  vs.  Millstone  Hill  1000  |xs  ISR  Te  Percent  Difference  Case  Summary 
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jendix  B;  DMSP  vs.  Millstone  Hill  TSR  T.  Comparison  VaIups 


Table  22.  DMSP  vs.  Millstone  ISR  Te  Comparison  Values 


Case  # 


Space¬ 

craft 


Zenith  Avg 
Angle  Kp 


DMSP  Avg. 
Te(K) 


2921.0 


2809.7 


4988.4 


4522.9 


4716.7 


4356.7 


3211.5 


4071.7 


3856.4 


3412.8 


2960.5 


4146.2 


4223.8 


2903.1 


4404.2 


2702.6 


410ps 
Ext.  Te 


2207.91 


2957.98 


2600.31 


2890.01 


3088.92 


2905.32 


2530.55 


2225.94 


2037.91 


3252.33 


2919.14 


2457.66 


3065.20 


1985.81 


3095.07 


2773.87 


2219.99 


2888.39 


2551.80 


3158.72 


2428.68 


2349.75 


3264.78 


2548.46 


2742.37 


3205.74 


2511.42 


3283.46 


3446.84 


3010.32 


2305.00 


3105.89 


2801.96 


2325.14 


2161.87 


1000  ps  Avg.  2000  ps  Avg. 
Avg.  Te  Error  Avg.  T.  Error 


Notes:  *  N/A 

*  Based  upon  only  two  profiles 
^SN  ratio  «  0.1 

13  :  Conjunction  in  Dark  (Terminator  Zenith  Angle  =118°) 
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